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Abstract 



The problem of determining the effects of the surrounding plasma on nu- 
clear reaction rates in stars is formulated ab initio, using the techniques of 
quantum statistical mechanics. Subject to the condition that the nuclear 
reactions ensue only at very close approach of the fusing ions and the con- 
dition that the reaction be slow, the authors derive a result that expresses 
the complete effects of Coulomb barrier penetration and of the influence of 
the surrounding plasma in terms of matrix elements of well defined operators. 
The corrections do not separate into the product of initial state and final state 
effects. When the energy release in the reaction is much greater than thermal 
energies, the corrections reduce, as expected, to evaluation of the equilib- 
rium probability of one ion being very near to the position of another ion. 
We address the calculation of this probability in an approach that is based 
on perturbation theory in the couplings of the plasma particles to the two 
fusing particles, with the Coulomb force between the fusing particles treated 
nonperturbatively and interactions among the plasma particles treated in the 
one-loop approximation. We recapture standard screening effects, find a cor- 
rection term that depends on the quantum mechanical nature of the plasma, 
and put an upper bound on the magnitude of the further correction terms 
for the case of a weakly coupled plasma. We find that possible "dynamical 
screening" effects that have been discussed in the literature are absent. The 
form of our results suggests that an approach that relies on numerical cal- 
culations of the correlation functions in a classical Coulomb gas, followed by 
construction of an effective two body potential and a quantum barrier pene- 
tration calculation, will miss physics that is as important as the physics that 
it includes. 
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I. INTRODUCTION AND SUMMARY 



The theoretical determination of fusion rates in stars [For example, Clayton (1968)] 
divides into two steps: I. Finding the nuclear "S factor" for specific reactions, which, loosely 
speaking, characterizes what the rate of a specific nuclear reaction would be in the absence 
of the Coulomb repulsion between the fusing ions. II. Taking into account both the Coulomb 
repulsion between the fusing particles and the influence of the surrounding charged plasma 
[Reviewed in Ichimaru (1993)]. There is quite a clean division between these two steps under 
conditions such that the energy release in the fusion is very large compared with thermal 
energies, and such that the reaction rate is small. Both of these conditions are well satisfied 
in the case of nuclear reactions in the solar interior. In such cases, the result that we seek 
in step II is, to first approximation, the value at zero separation of the equilibrium density- 
density correlation function of the fusing ions.Q We shall say nearly nothing further about 
step I in this paper, and turn directly to the barrier penetration and plasma aspect of the 
problem. 

The history of this problem begins with the work of Gamow (1928), who calculated the 
Coulomb barrier effects in the limit in which the only Coulomb force taken into account is 
that between the two fusing bodies, which gives the answer for the limit of an extremely 
dilute plasma. The first important corrections to the Gamow rates are found in the results 
of Salpeter (1954) for the effects of Debye screening by the electrons and ions in the plasma. 

1 The fusion reaction itself, and the associated nuclear forces, change the behavior of wave functions 
at very short relative distances. But as long as the reaction rate is very small, either because of 
the Coulomb barrier penetration factor, or because the fusion interaction is weak, these changes 
are significant only over a very small volume, and the exterior plasma physics will have negligible 
effect on the space variation of wave function within the volume. The quantities that determine 
the influence of the plasma are then the correlation functions that would exist with the nuclear 
charges concentrated at points and the nuclear interactions turned off. 
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These corrections, which will serve as the benchmark for comparison with further corrections, 
are given, at temperature?] T = by 

T 5 = r exp{/?eie 2 K D } , (1.1) 

where To is the rate in the absence of the plasma, e\ , e 2 are the charges of the fusing nuclei, 
and Kd is the Debye wave number. This wave number squared is the sum of those for each 
species in the plasma, 

«i> = (1-2) 

Here we count each spin state of a particular kind of particle with a separate species index s. 
The Debye wave number is generally defined by the rate of exponential fall off of the charge- 
charge correlation function of the plasma. For a dilute plasma, the one-loop or random 
phase approximation gives, as is reviewed in Appendix A [Eq. (|A32 )], 



4 ims fJ^k t (1 . 3) 

where (n s )j3 is the average density of species s which has charge e s = Z s e, and /i s is the 
chemical potential of this species. The general form may be needed for the electron compo- 
nent in the plasma of a star in which the effects of Fermi statistics may be important, but 
for the ionic species, Boltzmann statistics suffices, in which case the Debye wave number 
assumes the familiar form 

k\ s = iire 2 J(n s )p. (1.4) 

The effects of the plasma on reaction rates in the sun, as estimated from the above 
formulae, are modest but important. For example, in the entire region from the center of 
the sun to the radius inside which one-half of the energy is generated, the screening effects 
give about a 5% increase in the rate of proton-proton fusion and a 20% increase in the rate 



2 We measure temperature in energy units so that the Boltzmann constant is unity, fejg = 1. 
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for the reactionQp+ 7 Be — > 7+ 8 B. In lower temperature main sequence stars, and in highly 
evolved stars, the corrections become even more important. 

It is worthwhile exhibiting the physical origin of the Salpeter correction shown in 
Eq. ( | . 1| ) , the leading correction in the dilute plasma case. A plasma shields a charge 
brought into it, with a point charge q producing the Debye potential 

0(r) = ^ e - KDr . (1.5) 

Thus, assembling a charge distribution p(r) in a plasma has an associated polarization energy 
- the energy to assemble the charge less the corresponding vacuum energy — given by 

"g-Kolr-r'l y 



e P =\j J{dv){dv')p{v) 
The limit for a point charge q yields 



r — r r — r 



p(r'). (1.6) 



\<?k d . (1.7) 



The Boltzmann factor associated with this polarization energy alters the number density 
relation for the ions from that of a free gas to read 

(n) p = A- 3 e^e-^ , (1.8) 

where A is the thermal wavelength defined by 



with M the mass of the particle. A nuclear reaction rate depends upon the probability to 
find two particles at the same spot, or more precisely, upon the average of the product of the 
number densities of the two particles, (ni(0)n 2 (0)). With the two nuclear particles having 
charges t\ and this is given by 



3 These estimates use the tables in Bahcall (1989), p. 90. The density, H/He ratio, and temperature 
all change considerably over this region, but the screening factor almost remains constant. 
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(ni(0)n 2 ((% = Ar 3 A 2 -V^ + ^)exp{-/?e (ei+e2) } . (1.10) 

Thus, in view of the single-particle number density relation ( |1.8| ) and the form ( |1.7| ) of the 
polarization energy, we find that 

(n 1 (0)n 2 {0))/3= {n 1 )p{n 2 ) l3 e^{l3eie 2 ^D} ■ (1-11) 

The exponential is the Salpeter correction in Eq. ( |1 . 1[ ) . 

The method that we have just described is in the spirit of our work in this paper which 
makes use of grand canonical ensembles that entail chemical potentials. To make comparison 
with other treatments, we should note that the same result is obtained if one replaces the 
Coulomb potential between the fusing particles with a Debye screened potential of the 
form of Eq. (|1.5|) . Since only distances that are short on the scale of the Debye length 
Hp enter into the quantum tunneling, only the short-distance correction 5V = —e\t 2 ^D 
to the Coulomb potential e\e 2 jr need be retained. This is equivalent to an energy shift 
E — > E + eie 2 K£> or a momentum alteration Sp = eie 2 Kr,m/p. Thus the Gamow tunneling 
factor exp{— 27reie 2 m/|j} is corrected by the factor 

exp L ie2 K D r 7re ^ 3 2m j | • (1.12) 



In the thermal average, as we review below [c.f. Eq. ( 3.15 )1, the relative momentum p of 



the fusing particles is replaced by the most probable momentum p, where 

p 3 = 2ire 1 e 2 m 2 /{3 . (1.13) 

This replacement produces the Salpeter correction in Eq. ( |1 . 1| ) . 

The "weak screening" domain is defined as that in which the exponent in Eq. ( |1 . 1| ) is 
small compared to unity. This is also the condition under which the Debye formula gives a 
good approximation to the purely classical plasma correlations. Clearly, one way of going 
beyond the Salpeter correction and toward the domain of the "strongly coupled plasma" is to 
use a correlation function that is better than the one provided by the Debye formula. There 
is an extensive literature devoted to the computation of the classical correlation function for 



a plasma, and the determination of fusion rates therefrom [Salpeter and Van Home (1969); 
Dewitt, Graboske and Cooper (1973); Graboske et al. (1973); Jancovici (1977)]. In these 
approaches one translates a classical correlation function, usually numerically determined, 
into a modification of the potential between the fusing particles, and then calculates the 
quantum mechanical barrier penetration factor, using this potential, in order to obtain the 
fusion rate from a Gamow factor appropriate to this potential. We call this the "basically 
classical" approach. 

In this approach, there are conceptual problems raised by the division of the problem 
into a quantum mechanical and a classical part. The literature lacks any development that 
begins with a correct general expression for the rate, shows how a division into classical 
correlation plus quantum tunneling can be made as an approximation, and gives a system 
for finding the corrections to this approximation.^ 

We shall develop a general formulation of the problem which is based on only one ap- 
proximation — that the fusing nuclei in the plasma can themselves be described by Maxwell- 
Boltzmann statistics. In practice, this is an excellent approximation as these particles are 
seldom so dense as to require quantum statistics. No approximation is required for the 
remainder of the plasma. This basic approximation enables us to disentangle the internal 
Coulomb corrections in the fusion process from the Coulomb interactions of the fusing nuclei 

4 The work of Alastuey and Jancovici (1978) does not completely meet these criteria for at least 
two technical reasons: 1) Dividing the system into subsystems of fusing particles and surrounding 
plasma, the authors treat the plasma completely classically, thereby missing the correction term 
that turns out to be numerically most significant for the weakly coupled case. 2) Their approxi- 
mation scheme in which a numerically determined classical correlation function for the plasma can 
be applied to the quantum tunneling calculation for the fusing particles involves a zeroth order in 
which the fusing particles are frozen in an "average position" in space. As a result the effects of 
center of mass motion are irretrievably lost. 
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with all the other particles in the plasma (including other nuclei that could fuse). We then 
employ this general formulation to obtain a straightforward approach to the problem, based 
on perturbation theory. Considered as a factor, the Salpeter correction differs from unity 
by a term of leading order e 3 and therefore should be calculable in a perturbation expan- 
sion. In describing our approach it will be convenient to identify two particular ions as the 
"fusing particles" and all of the rest of the ions and electrons in the surrounding medium 
as the "plasma particles". As we shall see explicitly later, the Salpeter correction can be 
derived as the leading correction from a perturbative treatment in which the Coulomb force 
between the fusing particles is treated exactly, the coupling between the fusing particles and 
the plasma particles is treated perturbatively, and the plasma itself is treated in one-loop 
order. The formulation does not divide the development into a quantum mechanical and a 
classical part. 

We shall develop this approach in some detail, and address the corrections beyond the 
leading Salpeter term. We are still limited in application to a weak screening domain of 
density and temperature, but we are able to compare the correction terms with those that 
can be produced in the basically classical approach. If we had found rough agreement, then 
we could have argued that the basically classical approach is a reasonable way to approach 
the regions of stronger coupling. However we find the opposite, that the leading correction 
terms that emerge from the complete perturbative approach have a nature that is essentially 
different from those produced by the basically classical approach. Our results thus point to 
the need for a computational framework that does not divide into a classical and quantum 
part. 

The new quantitative results applicable to the weak screening domain are: 
1) It is frequently stated that the Salpeter result, or static screening, should be a good 
approximation only when the velocities of the plasma particles are greater than the veloc- 
ities of the fusing particles, so that the latter see an average potential rather than a time 
dependent one [e.g. Johnson et al. (1992), Shaviv & Shaviv (1996).] This condition is rather 
badly violated for the ionic component of the plasma. But we find, in the leading order 
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of e 3 , no "dynamical" modifications of the Salpeter result in the ionic component of the 
screening, in fact no dependence on the ionic masses at all. This conflicts with the results of 
Carrero, Schafer, and Koonin (1988), who found such terms. These "dynamical" corrections 
to screening, if present, would produce moderate decreases in predicted fusion rates under 
solar conditions. In Appendix D we show how such terms, which never arise in our approach, 
are produced in a real time treatment. We also find the terms that cancel them. 

2) There is one correction at the e 4 level that is numerically significant under solar 
conditions. The correction comes from a quantum mechanical term in the electronic part of 
the plasma response function, and is inaccessible to approaches based on a classical treatment 
of the plasma. 

3) A bound is put on all remaining corrections at the e 4 level. This bound limits the 
contributions of such terms to around the 0.1% level, under solar conditions. 

The sum of our new corrections turns out to be of marginal importance to the various 
solar neutrino puzzles [Bahcall (1989), (1995)]. Thus we conclude, as have others, that the 
resolutions to these puzzles will not be found in plasma physics. 



takes place over a very short distance in comparison with the particle separations in a 
plasma. Hence, in the non-relativistic, second-quantized formalism which we employ, it can 
be described by an effective local Hamiltonian density 



II. REACTION RATE THEORY 



A. General Framework 



A nuclear reaction, which we schematically indicate by 



1 + 2^3 + 4, 



(2.1) 



H(r,t) =/C(r,t) + /C t (r,t), 




in which 
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/C(r, t) = g e~ iQt ^(r, t) 4(r, t) ^ 2 (r, t) ^(r, i) , 



(2.3) 



where Q is the energy release of the reaction. The number density operator of a produced 
particle, say the particle labeled by 4, is given by 



n(0) = 4(0)^4(0) 



(2.4) 



Neglecting the kinematical particle flow described by the divergence of the corresponding 
particle flux vector which vanishes in the ensemble average, the time rate of change of this 
density is given by 



n(0) = -ij (dr) [n(0), H(r, 0)] 
= -i{/C(0) -/C f (0)} . 



(2.5) 



These operators create and destroy particles, and hence their expectation vanishes in the 
unperturbed plasma ensemble which is diagonal in particle number. The reaction rate T 
appears in the additional linear response to the action of the perturbation 7i, 

r0 



r = -i 



dtJ(dr)([h(0),H(r,t)])p 
= - f dt y'(rfr)([/C(0)-/C t (0),/C(r,t) + /C t (r,t) 



(2.6) 



where (• • -)p denotes the grand canonical thermal average of the background plasma. 
Since the operator /C does not conserve particle number, 



£(0),/C(r,f) 







ICt(0),/Ct(r,*) 



(2.7) 



Moreover, in virtue of the space-time translational invariance of the ensemble and the anti- 
symmetry of the commutator, 



/C( ),/Ct(-r, -t) 



P 



Thus, changing the sign of the integration variables in this term, (r, t) 
combine to produce an integral over all times, 



(2.8) 



[— r, — t), terms 
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r = -J_™dtJ(dr)([)C(0),tf(r,t)]) p . (2.9) 

The two terms which comprise the commutator which appear here correspond to the de- 
struction of particle 4 via 3 + 4-^1 + 2 as well as its creation by 1 + 2 — > 3 + 4. We shall 
assume that initially there are no particles of type 3 or 4 present in the plasma — or at 
least that they are extremely dilute — so that the operator K.'(r,t) acting to the right on 
an initial state vanishes. Hence, in this case, the only case that we shall consider, 

+0 °dtJ(dr) (£t(r,t)/C(0)) 



— oo 



g 2 f + °°dt /(dr)e^(v4(r,*)V4(M)V^ , 



oo 



(2.10) 

The result ( |2.10[ ) is equivalent to the 'golden rule' rate formula 



r = E^GOI ( F = 3 ' 4 i /|£(0)|J = 1, 2; i)\ 2 (2k)%E f - E I - Q)6(P F - Pj) . (2.11) 

Here \I = 1,2; i) describes the complete initial state of the system — the initial particles that 
take part in the reaction are labeled by 1, 2 and the background plasma by i. The initial 
particles 1, 2 are interacting with this background plasma, and so their separate energy- 
momentum is not conserved. The initial states, however, can be chosen to be eigenstates 
of the total energy-momentum Ej , P/, and this we have done. The probability distribution 
for the initial states is denoted by P(I)- Similarly, the final state (F = 3,4; f\ contains the 
produced particles 3, 4 and describes the final plasma / with the total energy-momentum 
of this state denoted by E F ,P F . The fusing particle coordinates and plasma coordinates 
are interlocked in a complex way in the states, making the golden rule formula somewhat 
ill adapted to proceeding further. In applications in which reaction energy release Q is very 
large compared with thermal energies, or in a case in which the final particles 3 and 4 are 
both neutral, it is possible to express the rate in terms of a thermal correlation function of 
particles 1,2 in the initial state, but we can proceed quite some way from ( 2.101 ) without 
making these restrictions. 
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B. Free Gas 



We pause now in our general development to illustrate the nature of the result (|2.10|) 
which we have just derived by considering the limit in which the plasma is replaced by a 
free gas. For a free gas, the ensemble expectation value in Eq. ( 2.10 ) factorizes. Since the 
particles 3 and 4 are produced out of the vacuum, one factor is 







^ 3 (r,t)^(r,t)Vi(0)Vl(0) 







(jgg) f ( rf P 4 ) J( P3 +P4)-r-i(i?3+£4)t 

(27r) 3 7 (2tt) 3 



(2.12) 



where 



The other factor is 
where 



2Af 3 



(dpi 



E A 



13 J (2vr) 3 



2M 4 



(dp 2 
(2k) 



(2.13) 



r n 2 (p 2 )e 



-i(pi+p 2 )-r+i(E 1 +E 2 )t 



(2.14) 



and 



1 



(2-15) 



^2(P2j 



(2-16) 



e P(E 2 -H2) ZSZ 1 

are the Bose or Fermi distributions of the initial particles with chemical potentials fii and 
/i 2 and energies 



(2.17) 



2Mi 2M 2 

With these results in hand, it easy to see that in the free gas limit, the rate formula ([2.1 0|) 
yieldsQ 



5 Here we have neglected possible spin variables. Spin degeneracy is readily accounted for by 
including spin-weight factors in the density distributions ni(pi) and n 2 (p 2 ). 
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f (dp*) 


f(dp 3 ) f 


1 (2tt) 3 J 


(2tt) 3 J 



(2tt) 3 zv ^ y J (2tt) 3 lv ^ iy 
(27r) 3 5(p 4 + p 3 - p 2 - Pl )(27r)5(E 4 + E 3 - E 2 - E l - Q) . (2.18) 

We use this free gas result to make several points. First of all, we note that it reproduces 
the exact nuclear rate 

if the coupling g 2 is taken inside the integration and replaced by the squared nuclear matrix 
element |T| 2 . In practice, the effective Hamiltonian method produces the correct nuclear rate 
if the coupling g 2 is replaced by an appropriate average value of the squared matrix element 
\T\ 2 . Note that, as yet, the rate formula does not account for the Coulomb interactions 
between the reacting particles. Finally, we note that for most applications such as the 
plasma in ordinary stars, the initial particles 1 and 2 are in a dilute gas. In this case, the 
momentum-space densities may be replaced by the Maxwell-Boltzmann forms, 

n( Pl ) = e -/^i-w) ? n ( p2 ) = e -/3(E 2 -M 2 ) _ (2.20) 

Expanding the Bose or Fermi denominators, we see that the correction to this leading 
approximation involves the very small factors exp{/3/ii} and exp{f3fi 2 }- Since in the dilute 
gas limit 

{nh = I^$ = A ~ V " ' (2 ' 21) 

we see that the corrections to the dilute limit are of order 



e 



(n^A 3 , (2.22) 



which, for example, at the core of the sun is roughly of order 10~ 6 . The corrections to the 
Maxwell-Boltzmann limit involve exchange effects of quantum identical particles. We turn 
now to exploit the dilute nature of the initial reacting particles so that they may be treated 
with classical statistics. 
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C. Plasma 



We shall work out the nuclear reaction rate accounting exactly for the Coulomb interac- 
tions between the reacting particles but taking advantage of their assumed dilute character 
so that we may neglect Bose or Fermi exchange effects. The other particles in the plasma, 
however, — such as electrons — need not be dilute, and they will not be restricted to obey 
classical statistics. 

We shall assume that the produced particles, which are labeled by 3 and 4, are sparsely 
distributed in the plasma so that their initial state which appears in the rate formula (ETTU) 



thermal expectation value (• • -)p can be represented by the vacuum state. Thus the evalua- 
tion of the rate involves the calculation of the vacuum expectation value in the subspace of 
the produced particles 

W(r 3 ,r 4 ,t) = (03,0 4 |^3(r3,t)^4(r4,t)^(0)^(0)|0 3 ,0 4 ) (2.23) 

in the limit r 3 = r = r 4 . The Heisenberg equation of motion for the field operators gives 

Ao 3 , 4 | V 3 (r 3 , t) </> 4 (r 4 , t) = (--3- - -S- + 6364 ) (0 3 , 4 |^ 3 (r 3 , t) ^ 4 (r 4 , t) 
at { 2M 3 2M4 |r 4 — r 3 | J 

+ (0 3 , 4 |V 3 (r 3 , t) Mu, t) [e 3 0(r 3 , t) + e 4 0(r 4 , t)] . (2.24) 

Here the operator <p(r, t) is the potential produced by all the other particles in the plasma. 
It is determined by the operator p(r',t) which measures the charge density produced by all 
these particles, 

0(r,t)= [(dr')—!—p(T',t). (2.25) 
J |r — r I 

In the absence of this potential, the bra vector in Eq. Q2.24j ) defined by the field operators 
acting to the left on (0 3 , 4 | obeys 

d 

^(r 3 , r 4 , t\ = (r 3 , r 4 , t\H 3A , (2.26) 
where, with ordinary quantum mechanical notation, 
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H 



3,4 



P§ + P4 



e 3 e 4 



2M 3 2M 4 |r 3 -r 4 | 



^2.27) 



Thus, 



r 4 , *l = ( r 3 



^2.28) 



The addition of the potential terms can be included by passing to an interaction picture 
with the time evolution operator 



Wit) = (expj-i f dH [e 3 <P(r 3 (t'),t')+e 4 <P(n(t'),t')]j 



a time-ordered exponential in the real time t with 



r 3or4 (t) — e lH3At r 3or4 e 



-iHzAt 



Thus 



(0 3 ,0 4 |V 3 (r 3 ,t)V'4(r 4 ,t) = (r 3 ,r A ,t\W(t) 



and we conclude that 



^2.29) 



(2.30) 



(2.31) 



W(r 3 ,r 4 ,t) = (r 3 ,r 4 ,t|^(t)|r 3 = = r 4 >, 



(2.32) 



since this obeys the equation of motion (|2.24|) and the correct boundary conditions when 
t — > 0. The rate formula ( |2.10|) is now expressed as 



dt I (dr) e lQt (4(r, t) ^(r, t) W(r, r, *) ^ 2 (0) ^i(O) 



(2.33) 



It should be noted that W(r, r, t) is an operator in the plasma state space. 

The thermal expectation value which appears in Eq. ( [2.33D is defined by a sum over all 
states (except for the final particle states which have been accounted for), a trace weighted 
by the density operator 



V = jex V {-(3 [#-X>«iV«]} 



^2.34) 
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where N a are the number operators of the particles with chemical potential //„, and Z is the 
partition function 

Z = Trexp{-[3[H -^aNa}} . (2.35) 

The density operator V defines the grand canonical ensemble that we employ. Since an 
exact treatment of the Coulomb interactions of the particular reacting particles requires 
an emphasis on the individual particle aspects, we are forced to employ the coordinate 
representation for these particles. Thus we use the coordinate representation for these 
particles in the trace that defines the average in the grand canonical ensemble and write, in 
a necessarily rather schematic way, 

(^(r, t) 4(r, t) W(r, r, t) ^(0) ^(0))^ = £ ^ J ■■ ■ 

(r[... ; r' 2 ... ; • • • \V ^(r, t) 4(r,t) W(r, r, t) ^(0) MVK • • • ; r' 2 •••;•••) . 

(2.36) 

The sum is a sum over all particle numbers N\ , N 2 (which are now just numbers, not 
operators). The coordinates of all the particles that appear in the state vectors are to be 
integrated over all space. The abbreviated notation used to label these particles in a state 
vector such as 

|r;...;r' 2 . ••;...> 

is as follows. The first set • • • denotes the coordinates of the Ni particles of type 1 ; 
the second set r' 2 • • • denotes the coordinates of the N 2 particles of type 2 ; the remaining 
• • • stand for all the other particles. The sum defining the trace for these other particles 
may be done in any representation. In view of the symmetry in the integration over the 
initial reacting particles coordinates, the action of the destruction and creation operators is, 
effectively within this integration, 

MO) MWi • • • ; r 2 •••;•••)= I ; )N 1 N 2 5(r[)5(r> 2 ) . (2.37) 

We now make the restriction that the initial reacting particles 1,2 are sufficiently dilute 
so that Boltzmann or classical statistics can be employed. In this approximation, exchange 
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effects are neglected, and the particles that were destroyed by the field operators in Eq. ([2.371) 
are identified with the particles destroyed when the remaining operators ip\(r, £) and ip 2 (r, t) 
act to the left. Accordingly, the ensemble average which appears in the rate formula (|2.33[) 
may be expressed as 

(4(r, t) 4(r, t) W(r, r, £) V*(0) ^(0))^ = £ ^ \ ^ ^ 1} , / • • • 

(•■■;■■•; ••■|(r' 1 = = r / 2 |^ 1 t (r,t)^(r,t)W(r,r,t)|0 1 ,0 2 )|--- ; •■■ ; • ■ ■) 
= Tr^ = = r' 2 \v4(r,t)4(r,t) |0 X , 2 )W(r, r, £) . (2.38) 

The action of the creation operators on the vacuum state |0i, O2) which appears here is 
essentially the adjoint of that appearing in Eq. ( |2.31| ), and so we have 

^(r, £) 4(r, £) |0 l5 2 > = V(t)\r[ = r = r' 2 ,t). (2.39) 

Here the time dependence of the state is governed by the quantum-mechanical Hamiltonian 

Hl ' 2 ~2M 1 + 2M 2 + ' (2 ' 40) 

with 

|ri = r = r 2 , t) = e^ 2 *|ri = r = 4) , (2-41) 

and 

V{t) = (expj* J'dt' [e 1 0(r 1 (t'),f) + e 2 0(r 2 (£'), £')]})_ (2-42) 

is an anti-time-ordered exponential in the real time £ with 

rio r2 (£) = e i ^r lor2 e- i ^. (2.43) 

To deal with the action of the statistical density operator V, we write the total Hamiltonian 
for the whole system which appears in V as 

H = H 1>2 + ei 0(n) + e 2 0(r 2 ) + H R , (2.44) 
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where 0(r) is the potential that the remaining particles in the plasma produce at the point 
r, and is the Hamiltonian for the remaining particles in the plasma. In this way, we can 
make use of the interaction picture in imaginary time in writing 

exp{-(3H} = exp{-/3(# 1>2 + H R )} U{f3) , (2.45) 

in which 

U{(3) = (exp j- ^dr[e 1 0(r 1 (r),r) + e 2 0(r 2 (r),r)]|j , (2.46) 

is a time-ordered exponential in the imaginary time r with 

ri or2 (r) = e H ^ T r lov2 e- H ^ T . (2.47) 

The explicit r dependence in the potential terms appears because the charge density of the 
background plasma moves in imaginary time under the action of H^. That is, 

0(r,r)= [{dv')—^— P {v',T), (2.48) 
J |r — r | 

where 

p(r', t) = e HRT p(r')e- HRT . (2.49) 



Using Eq. ( [2.45|) , we explicitly separate out the initial particles from the density operator 
to obtain 

<p = e P(^+^)p e -P H ^u{(3) , (2.50) 
where V is the density operator for the remaining particles. Hence with 

(r; = = r 2 |e"^- 2 = (r; = = r 2 , -i(3\ , (2.51) 



we have 



(ri = = r' 2 \v4(r,t) 4(r,t) |0 1} 2 ) = 

expiPfa + fi 2 )}P (ri = = r' 2 , -i(3\U {(5)V {t)\v' x = r = r' 2 , t) . (2.52) 



Inserting this result in Eq. (|2.38|) produces a plasma thermal expectation value, and so 
Eq. (F33D yields 

/ + OO r 
dt / (dr) e iQt 
-CO J 

'(ri = = r' 2 , -«/3|l7(/3)V(t)|ri =r = r' 2) t) 



= r = ri,t|W(tK = = 14))^. (2.53) 

This expression for the reaction rate is the major result of this paper. It separates out the 
plasma interactions with the particles undergoing the nuclear reaction from their internal 
Coulomb interactions. The plasma interactions have two different effects. One is to alter the 
thermal distribution describing the density of states of the initial two particles that enter 
into the reaction. This alteration is described by the operator U {(3). The other type of effect 
of the background plasma is to alter the dynamics of the motion of the particles undergoing 
the nuclear reaction. This is described by the operators V(t) and W(t). 

Since our notation is rather condensed, it is worthwhile describing the meaning of this 
result in a little more detail. Spelling out the time dependence of the states expresses 

dt / (dr)e iQt 

-oo J 

'(r[ = = r' 2 |e-^ 2 ?7(/3)F(t)e i ' ffl ' 2t |r / 1 = r = r' 2 ) 

(r' 3 = r = r' 4 \e- lH ^W(t)\ r ' 3 = = . (2.54) 

Consider, for example, the second matrix element that appears within the thermal expec- 



tation value. Writing out the explicit expression ( |2.29| ) for the operator W(t)), this is the 
matrix element 

(r 3 = r = <,tjW(t)|r 3 = = r 4 > = 

(r' 3 = r = r 4 , t\ (exp j-i J* dt' [e 3 0(r 3 (t'), t) + e 4 0(r 4 (t'), *')]}) W 3 = = r' 4 ) . (2.55) 

The potential 0(r 3 (t'),t') is a field operator in the space of states that describe the thermal 
ensemble denoted by the overall expectation value (• • •) j g. This field operator is evaluated 
at the position r 3 (t') which itself is an operator in the quantum- mechanical two-particle 
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states of the final produced particles. The expectation value of Eq. (|2.55|) , which is in the 



two-particle space of the produced particles, yields a functional of the field operators 4> that 
depends parametrically upon the coordinate r and the time t. A similar functional of the 
field operator <p is produced by the other, initial particle, matrix element in Eq. fl2.53| ). 



The reaction rate V is then obtained by the space-time integral of the plasma, thermal 
expectation value of these functionals as shown in Eq. ( |2.53| ). 

In solar fusion processes, the energy release Q of the nuclear reaction is much larger than 
the temperature T. In such cases, the integrand in the time integration in Eq. Q2.53Q is 



dominated by the rapidly varying phase factor exp{iQt} and the time dependence of the 
states (ig = r = r' 4 ,t\. Including only this time variation, and setting t = elsewhere 
in Eq. ( |2.53|) , the time integration produces an energy- conserving S function that sets the 
energy of the produced particles equal to the energy release Q. The time region that is 
important is of the order 1/Q, and the leading terms in the large Q limit are obtained by 
placing t = in V(t) and W(t), which are thus replaced by unity, and by placing t = in 
the state \r[ = r = r' 2 , t) as well, with corrections^ to this approximation roughly of order 
T/Q. We shall see this explicitly in later examples. Thus in the highly exothermic limit, 
the plasma dependence in the general rate formula ( |2.53| ) enters only in the quantity 

y i = r = r' 2 ,^\Umr' 1 = = r' 2 )) r 

(r' 1 = = r^|exp{-^ 1 , 2 } 
^exp|-| o /3 rfr[e 1 0(r 1 (r),r)+e 2 0(r 2 (r),r)]j^ J \r[ = r = r' 2 > . (2.56) 

The spatial coordinate integration over r in Eq. ( |2.53| ) produces momentum 5 functions that 
enforce the conservation of the total momentum in the nuclear reaction. As we shall also 



6 In the cases that we are considering, the thermal excitations in the plasma provide energy 
transfers of order T which is much less than Q. However, in the interiors of stars that are more 
highly evolved than the sun, less exothermic and even endothermic reactions become important, 
and then the full complexity of the time dependence in Eq. ( |2.53 ) must be taken into account. 
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describe later, the effect of this momentum conservation is to replace the ion masses in the 
plasma that have Coulomb interactions with the fusing nuclei by proper reduced masses that 
include the total mass of these nuclei. This takes into account the center-of-mass motion of 
the initial nuclear system through the plasma. 



III. RATE CALCULATIONS 



A. Dilute Limit 



We first turn to illustrate the character of our basic result (|2.53|) with the limit where 
the plasma is very dilute. In this limit, the Coulomb interactions of the reacting particles 
with the background plasma may be neglected, and thus the operators U(/3), V(t), and W(t) 
may be replaced by the identity operator. Thus in this dilute limit, the thermal expectation 
value is the trivial expectation value of the identity operator, and our result becomes simply 



/+oo /• 
dt / (dv)e iQt 
-oo J 



(ri = = r' 2 , -i(3\r[ = r = r' 2 , t) (r' 3 = r = r 4 , t\r' 3 = = r 4 ) 



(3.1) 



To evaluate the transformation functions that appear here, we introduce complete sets of 
intermediate states which are eigenstates of the total initial or final momentum, P or P', and 
of the initial or final relative momentum, p or p'. The relative momenta are the asymptotic 
momenta of the interacting scattering states. Thus 

(dP) 



and 







L 2> 



-z/3|r; 



4*> 



(2tt) 3 

(dp) 

(2tt) 3 



exp 
exp 



■£(/J-*)}hs£>(o)' 



(r 3 = r = r 4 , *|r 8 = = r 4 ) = J exp -z — t + zP • r 



(2tt) 3 
(dp') 



(27T) 



exp < —i 



P 



2m' 



^ } (0) 



(3.2) 



(3.3) 
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Here 



M .U| • M 2 .1/;; • M 4 

is the common total mass of the initial and final states, and 



111 111 



(3.4) 



(3.5) 



are the reduced masses of the initial and final states in the reaction. The squared amplitudes 
|-0 p (O) | 2 are the squared relative-motion Coulomb wave functions at the origin which have 
the form 

-2lT7] 



(3.6) 



with 



eie 2 m 


, e 3 e 4 m' 


V= | | , 


^ = i ,\ 


Ipi 


P 



(3.7) 

IPI IPI 
for the initial and final states. 

The space-time integration in Eq. ( |3.1| ) now produces energy- momentum conserving S 
functions, and we have 



(dP) \ P 2 
r c = # 2 exp{/?(/ii + /i 2 )} / 77^^ exp { -(3 



(27T) 



(2tt) 3 7 (2tt) 3 
Since Mm = M X M 2 . 



P 2 P 2 



(3.8) 



A X 3 A 2 3 



Mi \ 3/2 / M 2 \ 3/2 



„ , \ 3/2 / \ 3/2 



\2nf3j \2n(3 J \2-nf3) \2n i , 
and remembering the formula fl2.22|) for the number density, we see that 



(3.9) 



Tc = (ni)p(n 2 ) p g 2 



(dp') 



P 



3/2 



27rm 



J (dp) exp {-/^K^O) 



(27T) 



p' 



p^ 



2m' 2m 



(3.10) 
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This is of the form 



rc = (ni)p(n 2 )p / (dav)/3, 



(3.11) 



where now the differential cross section da for the reaction includes the Coulomb corrections 
for the initial and final states. 

Before turning to the plasma corrections, let us examine the scale — the size — of the 
variables that enter into the reaction rate calculation. Clearly, the common center-of-mass 
momentum for the initial and final states is governed by the Boltzmann factor. Thus, in 
order-of-magnitude , 



The corresponding extent of the motion of the center of mass in imaginary time is of order 



where A is the thermal wavelength of the center-of-mass motion. The magnitude of the initial 
relative momentum is determined by the balance of two factors, the Boltzmann factor which 
rapidly decreases as this momentum increases and the Coulomb barrier factor whose effect 
weakens as the relative momentum increases. In the usual applications, 2ttt] = 2neie2fn/p 
is a moderately large number in the momentum region of relevance, and the value (|3.6|) of 
the Coulomb wave function at the origin is dominated by the leading exponential so that the 
distribution of initial relative momentum is peaked at the minimum value of the exponent 
in the factor 




(3.12) 




(3.13) 




(3.14) 



which is the point p = p given by 



p A = 27re 1 e 2 m 2 /P . 



(3.15) 



At this point, the Coulomb parameter has the value 
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and the exponential displayed in Eq. ( |3.14j ) becomes exp{— 37r?y}. Note that 



A n = J^ (3.17) 
V m 



defines a thermal wavelength, while 

a 12 = (3.18) 

eie 2 m 

defines the Bohr radius of the initial reacting system. In terms of these parameters, 



(27T) 



2 



P=jri~- ( 3 - 19 ) 
a 12 a 12 

The corresponding kinetic energy p 2 /2m is, in general, much less than the energy release Q 
of the nuclear reaction. In these cases, the relative momentum of the produced particles is, 
in view of the energy-conserving 6 function in the rate ( |3.1U ), of order 



p' ~ pm'Q , (3.20) 

which is much larger than the typical initial relative momentum p. 

The value fl3.6|) of the Coulomb wave function at the origin is given, to leading order, by 
a simple tunneling process. The classical turning point r = r max which enters here, at the 
typical initial relative momentum p, is determined by 

f-H = 0, (3.21) 

or 

2eie 2 ^ 2fj 

'"max = ^5 = ~^ ■ (3.22) 

This sets the scale of the quantum motion in imaginary time for r(r). Indeed, the tunneling 
amplitude for the initial two-particle state in a thermal ensemble at temperature T = 1/(3 
can be obtained as the steepest descent approximation to the path integral representation 
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[Alastuey and Jancovici (1978)]. In this approximation, the exponential exp{— 3irf]} appears 
as the exponent of the classical action computed for a classical solution in imaginary time 
that has an orbit extending from r = to r = r max back to r = 0. This classical solution 
gives the leading approximation to the quantum motion. 

B. Plasma Corrections 

As we remarked in the discussion of the general rate formula (|2.53|) , for highly exothermic 
reactions we may neglect the effects of the operators V(t) and W(t) so that the only plasma 
dependence is contained in the operator U(0). The scale of spatial variations in the plasma is 
controlled by the Debye length which is much larger than the estimates of the size of the 
coordinates r lor2 (r) which was just presented. Hence, we may, in the first approximation, 
set r = in the potential operator 0(r(r),r). As we shall soon see, the imaginary time 
dependence in 0(r, r) may also be neglected in leading order so that we may approximate 

U(J3) = exp{-/3( ei + e 2 )0(O)} . (3.23) 

Assuming that the plasma ensemble may be described by uncorrelated, Gaussian statistics, 
we now have 

(UH3f))p = exp{/3 2 ( ei + e 2 ) 2 (l/2)(0(O)0(O)) /3 } . (3.24) 
Using the potential correlation function obtained in Appendix A, Eq. (|A39|) , 

(<Kr)<K* f )) P = t 1 - exp{-* D |r " r 'l}] , (3-25) 

we have 

(U(J3)) p = exp{(l/2)/5( ei + e 2 ) 2 K D } . (3.26) 

It should be remarked that the exponent that appears here is generally not a large number 
so that the Gaussian statistics assumption is not needed. As remarked in the Introduction, 
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factors combine to form the plasma corrected number densities, and as shown in Appendix 
B (assuming that the result exponentiates) 

(ni ol2 ) p = A^ r2 exp{^// lor2 }exp{(l/2)^e? or2 /«£,}. (3.27) 

Therefore, the leading plasma correction is given by 

T s = exp{/5eie 2 n D }T c , (3.28) 



where Yc is the Coulomb corrected rate (|3.10| ). And, as remarked in the Introduction, this 
is the Salpeter corrected rate. 

We turn now to examine the size of the terms that we have neglected. We consider 
first the "dynamical" correction. Since the important time region is the small time interval 
< \t\ < 1/Q, in which the reacting particles move little, in the first approximation we may 
neglect both the spatial and time coordinates in the potential operator 0(r, t). Hence, in 
leading order, the dynamical corrections involve 

V(t)W(t) = exp{it( ei + e 2 )0(O)} exp{-tt(e 3 + e 4 )0(O)} = 1 , (3.29) 

since charge is conserved in the reaction, e\ + e 2 = e 3 + e±. To obtain the first "dynamical 
correction" , we note that since short times are involved, the dominant term is that with the 
least number of time integrations. This term involves the interference of a correction from 
the "statistical" U((3) and one of the "dynamical" V(t) or W(t) factors, 

± z/3( ei + e 2 )e r jT* dt'(<f>(0) [0(r(f), t') - 0(0)])^ . (3.30) 

Here the ± is + for the V(t) factor of the initial particles 1,2 and — for the W(t) factor 
of the final particles 3, 4. The charge of one of the reacting particles is denoted by e r and 
r(t') is its position operator. This correction is examined in Appendix E which provides the 
estimate 
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where oo^ ~ n 2 D /f3M defines the plasma frequency of a typical ion in the plasma. We 
have reverted to ordinary units in displaying this factor with Planck's constant H to show 
explicitly that it is of a quantum mechanical nature. The factor in braces is generally a 
very small correction indeed to the basic Salpeter correction (/3eie2 «£>). For example, at 
the center of the sun, and for Q — 1 MeV, the factor in braces is about 1CT 6 . Although one 
might expect that the contribution of the electrons in the plasma with their larger plasma 
frequency would dominate the dynamical correction, as shown in Appendix E, the leading 
contribution of the electrons cancels when summed over all the reacting particles. 

We turn at last to estimate the error involved in the leading rate correction which comes 
from the "statistical" U(f3). (A related and simpler calculation of the single-particle number 
density is presented in Appendix B.) To make this estimate, we examine the initial particle 
expectation value of the dilute plasma limit 

((r' 1 = = r 2 ,^|^(/3)|r' 1 = r = r 2 ,t)) /3 
^(r; = = r 2 ,-^|jl + ijVdr' 

((W( ri (r), r) + e 2 0(r 2 (r), r)][e 1 0(r 1 (r / ), r') + e 2 0(r 2 (r'), r')]) + ) 

|r / 1 = r = r' 2 ,t). (3.32) 

The plasma expectation value that appears here may be written as a spatial Fourier trans- 
form 

((0(r, r)0(r', r')) + ) p = J ^^<^G,(k, r - r') . (3.33) 

To discuss the initial particle expectation value, we use the center-of-mass and relative 
variables, 

R= (M iri + M 2 r 2 )/M, P = Pl + p 2 , (3.34) 

and 

r = ri -r 2 , p = (M 2 p 1 -M 1 p 2 )/M, (3.35) 
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where, as before, M = M\ + M2 is the total mass of the reacting system. In terms of these 
variables this expectation value has the direct product form 







-2 • 



-i/3 



T [ = r = r' 2 , t 



R' = 



exp < —(3 



2M 







exp {— f3H r } ■ 



exp I %—t R' = r) exp {iH r t} f = O) , 



(3.36) 



in which 



P eie 2 
2m If I 



(3.37) 



is the Hamiltonian for the relative motion of the initial system with reduced mass m. 

We shall first work out the effects of the center-of-mass motion. In terms of the center- 
of-mass and relative coordinates, 



r! = R + M 2 r/M, r, 11 .U|i-/.\/. 



(3.38) 



Thus, using the Fourier transform ( |3.33| ) in evaluating the correction ( |3.32j ), we see that we 
need to calculate the time ordered product 



(exp{zk.[R(r)-R(r')]} 



The center-of-mass coordinate undergoes free particle motion in imaginary time, 

R( r ) = R-iPr/M. (3.39) 

Using the canonical commutation relations to order the exponential, 

exp{ik R + k • Pr/M} = exp{k 2 r/2M} exp{ik • R} exp{k • Pr/M} 

= exp{-k 2 r/2M} exp{k • Pr/M} exp{ik ■ R} , (3.40) 

we find that 

(exp{ik • [R(r) - R(r')]}) = exp{-k 2 |r - t'\/2M} exp{k • P(r - r')/M} . (3.41) 
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With this result in hand, we can compute explicitly the center-of-mass contribution to the 
initial particle expectation value needed for the correction fl3.32| ). Introducing a complete 
set of intermediate state total momentum eigenstates, we now have 

R' = 0| exp 1-/3^1 (exp {zk ■ [R(r) - R(r')]}) + exp ji^ij |r' = 



exp{-k 2 |r - t'\/2M} J i|^exp{-*P • r} 

(/3 - it)} exp{k • P(r - r')/M} . 



exp 



2M 



(3.42) 



The spatial (dr) integration which appears in the rate formula (|2.53|) produces a 5 function 
which identifies the initial and final center-of-mass momenta P and P'. Thus the time- 
dependent factor exp t\ in Eq. (|3.42|) cancels the factor exp j— t| which is contained 
in the final produced particle matrix element. Thus, effectively, we may delete all the time 
dependence associated with the center-of-mass motion, delete the r dependence, and remove 
the integration over the center-of-mass momentum P' of the final produced particles. With 
this understood, we may complete the square in the (dP) integration in Eq. Q3.42| ) and get, 
effectively, 



R' = 0| exp j -p— } ( exp {*k • [R(r) - R(r')]}) + exp j i—t 



2M 



R' 



M 



3/2 



in which 



C(k, 



T — T 



2ir(3l 



exp 



C(k,r-r'), 



(3.43) 



(3k 2 



2M 



T — T 



(3.44) 



where 



r — t 



T — T 



1 



r — t 



(3.45) 



P V P J 

In an attempt to restrain notational proliferation, we define the normalized expectation 
value of an operator X by 



X 







exp {— PH r } X exp {iH r t} 











exp {— PH r } exp {iH r t) 







(3.46) 
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We should note that the real time dependence that appears here — the t dependence - 
makes a negligible contribution, and one can set t = here. This time dependence only 
gives rise to energy alterations on the order of particle thermal energies, which are negligible 
compared to the energy release of the nuclear reaction. Using this notation and the results 
that we have thus far obtained, Eq. fl3.32| ) becomes, with the "effective" caveats noted above 



still in place, 

V 1 = = r' 2 , -iP\U(P)\v[ = r = r' 2 ,t) 



(r[ = = r' 2 , -iflrj = r = r' 2 ,t){l + i ^ drdr' J ^ G p(k,r - r')C(k,r - r' 



(| [e\ exp {tk ■ [f (r) - r(r')}M 2 /M} + e\ exp {-ik ■ [r(r) - y{t')]Mi/M} 
+2eie 2 exp jik 



M 2 Mi , 
— r r H r(r 

M M 



(3.47) 



To the order that we require, the potential correlation function Gp(k,r — r') is given by 
the "sum of ring graphs" approximation, as discussed in Appendix A. To get the leading 
terms for dilute plasmas, a simplification can be made by subtracting the value at k = 
of the function that multiplies Gp(k, r — r') in Eq. (|3.47| ), and adding it back as a separate 
term. In the difference term we can replace the "ring sum" for Gp(k, t — r') by the lowest 
order bubble. In this way we simplify the integration in a manner that avoids a potential 
infrared problem and captures all terms of order e 3 and e 4 in the final result for the rate. 
For the term in which the multiplying function is evaluated at k = 0, we retain the complete 
Gp(k, t — r') whose long-distance Debye screening removes the infrared divergence, and we 
encounter an expression evaluated in Eq. (|B16|) of Appendix B, 

f drdr 1 [ j^Gp(k,r-r')=PK D + 8 f drdr 1 f° dk \^°\k, r - r') . (3.48) 
Jo J (Ztt) 6 Jo Jo k dk 

Here the first term involving the Debye wave number k d , which is of order e not e 2 , isolates 

the long-distance contribution. The second term involves the first-order polarization function 

computed in Appendix A, whose dilute form reads 

n<°>(*, r-r')^ Y,/s (n.)f>C.(k, r - t>) , (3.49) 
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where 



C s (k,T-r') 



cxp 



(3k 2 
2M. 



fir - r') 



(3.50) 



Having removed the long-distance contribution, we may now make the replacement 



(3.51) 



in the remaining terms, which are of the form 



G p {k,T-T') C(k,r-r')e l - 1 



(3.52) 



The time ordering involves the operators f (r) and r(r'), each of whose components commute 
amongst themselves. Hence the operators in the time-ordered product may be treated like 
ordinary numbers in intermediate calculations. Thus we now encounter a sum of the form 



(3.53) 



of terms of the form 



(2tt)3 UV P l 2M s 



exp 



/(r-r')j 

!e f(r-r'))e^-l 



2M" 



(3.54) 



To evaluate this expression, we represent 



1 \ * /"00 



(3.55) 



and complete the square in the Gaussian wave number integrals to evaluate them and write 
the expression ( |3.54| ) as 

( v.2 



rco 

vAne 2 / sds 
Jo 



[s + a s ) 3/2 exp 



4(s + a s ) 



(s + b s )- 3 / 2 



(3.56) 



where 



JL 

2M. 



■f(r-r') 



(3.57) 



31 



with M s the reduced mass defined by 

1 

Ml 



1 1 

M + W, 



and 



2M. 



f{r-r'). 



(3.58) 



(3.59) 



The potential infrared divergence is reflected in the large s behavior of the integrand in 
Eq. ( |3.56| ) in which the separate terms do not give convergent integrals. The cancellation of 
the divergent pieces is exhibited by partial integration of the leading terms which reduces 
the expression in Eq. (|3.56|) to 



POO 

y/irel / ds 
Jo 



o (s + a s ) 3 / 2 



4a s + 



s + a s 



exp 



4(s + a s 



8eUhrbs 



(3.60) 



It is easy to check that the last term here involving Vb s gives contributions that precisely 
cancel the second term on the right-hand side of Eq. (|3.48|) . Omitting this term and making 
a simple change of variables presents the result as 



- 4e a y/7nT s 



1 du 
o \fu 



1 + 



u 



4a, 





ur 2 \ 


exp | 


4a J 



(3.61) 



As is discussed in Appendix B, the contribution of the electrons in the plasma should be 
separated out because their mass is so small relative to the nuclear reacting particles and 
the nuclei in the plasma. The relatively small electron mass means that effectively a s — > oo 
or that we may take r 2 = for their contribution in Eq. (|3.60| ). Thus, for the electron 
contribution, the reacting particles may be taken to be at the same spatial position, and the 
calculation of their effect is the same as that for the number density calculation presented 
in Appendix B, but for a particle of charge e\ + e2- We evaluate the electronic part of the 
second term on the right-hand-side of Eq. (|3.48|) using the results of this Appendix, which 
include the effects of Fermi-Dirac statistics [], to obtain 



7 Note that the Debye wave number which appears here should be computed with the quantum- 
mechanical form (1.3) for the electrons presented in Appendix A [( A32| )] if the Fermi statistics for 
these particles is important. 
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-2 > 



A , t) 



(r' 1 = = r 2 ,-i(3\r' 1 = r = r> 2 ,t) 



(3.62) 



in which 



( 



1 du 
u 



+e 2 2 



1 + 
1 + 



+ 2eie 2 



1 - u) 



- L f(r) - r(r')] 2 (M 2 /M) 2 
[f (r) - f (r')] 2 (M x /M) 2 



exp < — u- 



exp < — it 



[v{r )~v{r')f{M 2 /Mf 
4a s 

4a, 



4a s 



exp < — «• 



(r) + f f (r- 



4a e 



(3.63) 



As we shall soon see, the third term inside of the curly bracket in Eq. ( |3.62| ) is the most 
important correction, after the Salpeter correction, for the case of a weakly coupled plasma. 
It involves the quantum mechanics of the plasma in an essential way. This fact leads to a 
dilemma in considering how to approach the regime of lower temperature or higher density, 
a region in which the perturbation expansion is not valid. The standard approach to this 
domain is to use non-perturbative methods to treat the plasma classically, and then to use 
the extracted effective potentials in the tunneling calculation that determine the rates. But 
the fact that the leading perturbative correction to the Salpeter result requires a quantum 
treatment of the plasma argues that the quantum mechanics of the plasma will be central 
to any meaningful calculation of the strongly interacting case. 

The correction X due the ionic component of the plasma is quite small for a weakly 
coupled plasma, as a bound which we shall shortly derive proves. But before doing this, 
it is worth giving more detail to explain the remark of the previous paragraph that the 
character of this correction for the weakly coupled case brings into question the standard 
classical approach for dense plasmas. The matrix elements of the position operator f (r) are 
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controlled by the magnitude of the classical turning point r max discussed above [Eq. ( |3.22| )1. 



Thus the exponents in Eq. ( |3.63| ) involve the parametric ratio r^ ax /a s . In terms of the 
thermal wave length 

% = ¥r > ( 3 - 64 ) 

the denominator in this parametric ratio appears as 

a s = ^f(r,T>). (3.65) 

The thermal wave length X s is a quantum length — in ordinary units, it is of order Planck's 
constant h. At the center of a typical star like the sun, this quantum-mechanical parametric 
ratio is about unity. In such cases, no small parameter appears in the ionic correction, 
no further approximations may be performed, and a full quantum-mechanical evaluation is 
required. It must be emphasized that the relevant ratio r max / A s entails the time- dependent, 
quantum-mechanical aspect of the plasma. A static, classical treatment of the plasma for 
such corrections is unphysical and gives incorrect results. 

We now turn to bound the correction X and then evaluate this bound approximately. 
The bound is easily obtained from the remark that the quantum expectation values entailed 
in the correction X have path integral representations of the form 



F = 



cxp {-(3H r } (F(f (r), r(r')))Jr' = o) = J [dr}e~ s F(r(r), r(r')) , (3.66) 



in which the paths start at r(0) = and end at r(/3) = 0. The action functional S is 
the classical action continued to imaginary time. Since the action S is real, the measure 
[dr] exp{— S} is positive. Thus the expectation value ( |3.66| ) is bounded when the numerical 
function F(r(r), i"(r')) is bounded. Denoting an exponent which appears in Eq. ( |3.63| ) by 
ug, we note that the bounded combinations 

l-ug)e- U9 \<l (3.67) 

appears, giving the resulting parameter integral 
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1 du 



2. 



(3.68) 



Following the evaluation of Eq. (|B20| ) described in Appendix B, the imaginary time integrals 
for these terms involve 



7T 



drdr y/Aixa s = —f3 2 \ s . 

8 

Thus, bounding the remaining terms with 



(3.69) 
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(3.70) 



we obtain the bound 



1*1 <£/?«As{^(ei + e2) 2 A s + i f djiW 



A/3 2 Jo 



e ^ + e ^) |f(r) - f(r)l + 26162 
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(3.71) 



where 



"As 



47r/fe>„)$ 



(o) 



(3.72) 



is the contribution of the plasma species s to the squared Debye wave length. 

We should note that although the first term in the curly brackets in the bound (|3.71|) 
involves the quantum-mechanical nature of the plasma since it involves the quantum- 
mechanical wavelength X s , the second set of terms refer only to the classical aspect of the 
plasma, albeit as measured by quantum expectation values of the reacting particles. This 
second set of terms is, in fact, just the classical plasma limit of the ionic correction Q3.63| ). 
They are obtained from the result ( |3.63| ) if the limit A s — > is taken. This is the limit in 
which the plasma is treated classically, the formal h —>■ limit, or the limit in which the 
plasma particles are taken to have infinite mass (along with the total mass of the reacting 
pair). This limit is equivalent to taking a s — > 0. In the limit, the exponentials completely 
damp out the u integration except at the end point u = 0. Hence, the u integration may be 
extended to infinity. The terms which survive in the limit yield precisely this second set of 
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terms in Eq. ( |3.71 ). The connection of this limit with a computation which starts out with 



a classical plasma is spelled out in some detail in Appendix C. 

To estimate the size of the remaining quantum expectation values in Eq. ( |3.71| ), we 



note, as has been previously mentioned, that the basic quantum tunneling process is well 
described by the steepest descent approximation to the path integral. Thus we may estimate 
the size of these expectation values by replacing their quantum operators by the classical 
solution to the Coulomb problem in imaginary time which obeys the boundary conditions 
r(0) = = r(/?).0 It is easy to verify that this solution to 

d 2 r(r) eie 2 , , , , 

^ = -^r(r), (3.73) 

which gives the stationary path, has the parametric representation 

r = |^-sinO, (3.74) 

r ( r ) = ^r max k(l - cosf) , (3.75) 

where < £ < 2ir, r max is the classical turning point given in Eq. ( |3.22| ), and k is an 
arbitrary unit vector. The resulting integral 

/= f drdr \v(t)-x(t')\ 

= (£) ^max^^e-sinO^e'-sineOlcose-cose'l (3.76) 



8 The relative error to this classical evaluation is given by the first quantum correction which is of 
order Ar 2 (2/r max ) 2 , where Ar 2 is the average squared quantum fluctuation about the classical path. 
As an examination of the quadratic fluctuation correction to the path integral shows, this quantum 
fluctuation is of the same order as that for a harmonic oscillator with the same frequency, namely 
Ar 2 ~ l/(mu) = f3/{2irm). Using Eq's. fT22"D and ( p^ ), we find that Ar 2 (2/r max ) 2 ~ 1/fj, which 



is negligible when the Coulomb parameter rj is large. At the center of the sun, r\ ~ 2, and so our 
classical evaluation of the bound is a little rough. 
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may be computed by partial integrations which lead to 

d d 



dx dx' 



\x — X 



-28(x — x) . 



(3.77) 



Thus 



1 ^max 

3tt 2 



(3.78) 



The remaining integral is given by 



J = drdr' 
Jo 

M 2 + M 1 (3 1 
M 2n2 rT> 
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(3.79) 



Within this approximation we have 



\X\<J2f3K 2 Dy 



1 2 - 

— (ei + e 2 ) A s + r n 
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+ e2 Mj + 8 eie2 



(3.80) 



We return to our result (|3.62| ) and omit the correction due to the ionic component of the 
plasma. Using the Bohr radius of the hydrogen atom ao — 1/ e 2 m e and noting that the terms 
involving e\ and e 2 just provide the corrected number density - chemical relation which is 
described in Appendix B, we find that the plasma corrected rate is given by 

1 1 



1 + /3eie 2 



k d - 



c 



(3.81) 



2a e-^ + 1 

where Tq is the rate in the absence of plasma corrections given in Eq. ( |3.10|) . In writing 
this result we have neglected the quantum effects described by the correction X given in 
Eq. ( |3.63j ). For the example given by the center of the sun, r max / A ~ 1/2, and k d \ w 2x 10 -2 . 
Taking e\ = e 2 , we find that |X| < ^e\e 2 ^D x 1 x 10 -2 . That is, the quantum correction 
X is less than 1% of the basic Salpeter correction. At the center of the sun, the electron 
correction that has been kept in the result (|3.81 ) above, the second term in the square 
brackets, involves l/2ao ~ and exp{/5/i e } ~ 0.3. Thus it forms a correction of about 

10% relative to the basic Salpeter correction. Moreover, it should be kept in mind that with 
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exp{/5/i e } ps 0.3, the correct quantum expression (|1.3j ) for the electron contribution to the 
Debye wave number is reduced by about 10% relative to its classical value ( |1.4| ). 
Since the last, electronic correction which appears in Eq. ( |3.81| ), 

1 1 



K D 



(3.82) 



2a e~^» + 1. 

is the leading correction to the basic Salpeter result, we should describe its nature in a 
little more detail. In our formal counting of orders in the electric charge e, the Debye 
wave number kd is of order e, while the new correction is of order e 2 . The inverse Bohr 
radius 1/ao is proportional to one power of Planck's constant h, and so the new correction 
explicitly involves the quantum-mechanical character of the plasma. The final factor with 
the electron's chemical potential represents the effect of Fermi-Dirac statistics. When the 
electrons are sufficiently dilute (exp{/5/i e } <C 1) so that the Boltzmann limit can be used, 
the new correction may be written as 

1 



D,e /V e 



^3.83) 



Here k 2 d e is the electronic contribution to the squared Debye wave length, and A e is the 
electron's thermal wave length defined long ago in Eq. ( |1.9| ). The thermal wave length 
is again proportional to h, and we again see that the new correction is a first quantum 
correction. 



IV. CONCLUSIONS 



We have developed a concise expression [Eq. (|2.53|) 1 for the combined effects of the 



surrounding plasma and the Coulomb barriers on reaction rates in a completely ionized 
plasma. The expression is applicable to any case in which the actual reaction takes place 
at very small particle separations, and in which the reaction rate is slow, in a well defined 
sense. In general, there can be significant Coulomb and plasma effects in both initial and 
final states, and these effects are entangled in the results. For the cases of importance 
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in solar physics, in which the release of energy in the reaction is very large compared to 
thermal energies, the correction reduces to one which involves a factor [Eq. ( |2.56| )1 that is 
a generalization of the probability that one fusing ion is at the position of the other. This 
generalization takes account of the effects of the center-of-mass motion of the initial fusing 
nuclei through the plasma. 

We addressed the calculation of this quantity in an approach that is based on perturba- 
tion theory in the couplings of the plasma particles to the two (distinguished) particles that 
will undergo fusion. The Coulomb force between the fusing particles is treated nonpertur- 
batively through use of the action-minimizing imaginary time paths for tunneling problems; 
alternatively our expressions could be expanded in Coulomb wave functions. The interac- 
tions among the plasma particles were treated in one-loop order. Expressing the result for 
the rate as a correction factor times the result in the absence of the plasma, our calculation 
captures all terms in the correction factor of order e 3 and e 4 . 

The results consisted of a recapitulation of the standard Salpeter results for the effects 
of the electron and ionic plasma particles, plus a single new analytical term coming from 
the quantum nature of the electron plasma, plus a rigorous bound on the remaining terms 
of our expression, which turn out to make negligible contributions under solar conditions. 
We have two conclusions that are significant to calculations of solar processes: 

1) The quantum mechanical term in the dynamics of the electron plasma gives a reduction 
in the fusion rate of about 10 percent of the Salpeter enhancement term, under the conditions 
that prevail in the core of the sun. 

2) There is no "dynamical screening" modification of the Salpeter result, to leading order 
in the dimensional parameters. There is no dependence on the masses of the plasma ions 
of the terms that are appreciable under solar conditions. In consequence, the fusion rates 
are somewhat greater than those calculated by Carrero et al. (1988). In Appendix D we 



39 



demonstrate both how such terms arise in a real time approach and how they are cancelled.^ 
The nature of the remaining small correction terms in our results, and the formal struc- 
ture of the results as well, argue against an approach to the strongly coupled problem in 
which numerical results for the classical plasma are turned into effective two body poten- 
tials, followed by the quantum tunneling calculation to obtain rates. In the weakly coupled 
domain, however, the result of our systematic perturbation theory provides support for the 
way in which screening calculations have been used in the literature, subject to the modifi- 
cations mentioned above. In particular, the interparticle spacing does not have to be much 
smaller than the Debye length, even though this assumption goes into the derivation of the 
classical Debye potential.0 We reemphasize the fact that our treatment nowhere used this 
potential. 

The results of section 2 are general, and could be the first step in developing approxima- 
tions for the strongly coupled case. The perturbation development in section 3 is applicable 
only to the weak screening domain. In view of the smallness of the e 4 terms that we found, 
this fact seems to present no problems under solar conditions, but will probably not be 
useful for the interiors of much smaller main sequence stars or more highly evolved stars. 

9 As of the time of the writing of this paper it seems likely that the most critical calculation of 
rates in the sun, for the purpose of understanding the solar neutrino signal, will be the comparison 
of the p+ 7 Be — > 7+ 8 B rate with that of the rate of electron capture on 7 Be [Bahcall (1995), Hata 
and Langacker (1995)]. Modification (1) gives about a 2% decrease and modification (2) about 
a 2% increase for the p + Be reaction. Thus our results have little impact on the fusion side of 
the comparison. However, it would be interesting to evaluate the plasma effects on the electron 
capture process using the methods of this paper. 

10 We thus take issue with the remarks of Dar and Shaviv (1996), who question the applicability 
of the standard screening lore under solar conditions, based on the fact that this condition on the 
interparticle spacing is not satisfied. 
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APPENDIX A: THERMODYNAMIC FIELD THEORY REVIEW 

We review and collect here some results of multiparticle, non-relativistic field theory at 
finite temperature that are used in our work.Q 

The dependence in imaginary time of an arbitrary operator is given by 

X(r) =e HT Xe~ HT . (Al) 

The thermal average of two operators is defined by 

(A{t)B{t')) p = Z- 1 Tr exp {-/3 [h - E a ^A a ] } A(r) B(r') , (A2) 

where 

Z = Tr exp {-/?[#- E//a] } (A3) 

is the grand canonical partition function. Thus, if the operators commute with the particle 
number operators N a , there is the cyclic symmetry 

(A(t)B(t')) p =(B(t')A(t-(3)) p , (A4) 

which follows from the cyclic symmetry of the trace and the time dependence given in 
Eq. (|A1|) . Accordingly, the Green's function — the ordered product in imaginary time - 
of two operators, 



n The review here is close in spirit to discussion in Chapter 2 of Brown (1992), especially Problems 
3 and 4. 
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G(t-t0 = ((A(t)B(O) + ) , (A5) 
is periodic in imaginary time with period /3, and it may be expanded in a Fourier series, 

G(r-r') = -^E n ^ (T - r '^(^n), (A6) 

where 

u n = 2im/f3. (A7) 
Going over to real time (r — > it), with 

X(t) = e lHt Xe- iHt , (A8) 
we write the thermal average of the commutator as the Fourier integral 

( [A(t), 3(f)] >/3 = jT ^ c(u,) e -*-(^) . (A9) 



Using the cyclic symmetry (A3), a short calculation shows that, in the imaginary time 
interval < r < /3, 

Therefore, 

<7(iu; n ) = jyr(A(r)B(0)) p e^ T 

= r^^M_. (A11) 

J-oo 27r — zu; n 

To relate these results to the charge-density correlations that we need for our work, we 
note that first-order perturbation theory shows that the introduction of a small external 
potential 50 ex t(r, t) induces an average charge density given by 

(8p(r, t))p = J (dr') J dt'R{v -r',t- t') 50 ext (r', if) , (A12) 

where 

42 



R (r - r ' :t - t') = -i{ [p(r, t), p(r', t')} ) p 6(t - f 

(dk) r du 
(2tt) 3 J 2^ 



(A13) 



The similar Fourier transform of the commutator function itself, without the factor —i and 
without the step function 6{t — t'), defines a weight which we shall denote by c(k,u). Since 
this commutator function is odd under the interchange t <-> t', c(k, —uj) = —c(k,u) is odd 
in uj. Since complex conjugation of the expectation value the commutator of the Hermitian 
fields p is equivalent to this interchange, which is compensated by the complex conjugation of 
exp{— iuo{t — £')}, c(k, u)* = c(k, u) is a real function. The relationship between the Fourier 
transform of the retarded commutator function which appears in Eq. ( [A13j ) and the Fourier 
transform c(k, u) of the commutator function itself is given by the dispersion relation 

du' c(k,u') 



r(k, 00) 



(A14) 



-oo 27r uj — uj' — ie 

as one can easily prove by directly Fourier transforming over only positive time differences 
the Fourier representation of the commutator function. In view of the general relation (All) 
between the commutator and Green's functions, we see that the charge-density correlator 
in imaginary time, 



p(r,r)p(r',r') 



ft — 'r 



-iuj n {T-r') 



has a weight given by the analytic continuation 



(2tt) 3 



g{k,iu n )e 



ik-(r-r') 



(A15) 



g(k,iu n ) = -r(k,iu n ) . 



(A16) 



The induced charge - perturbing potential relation in Fourier space, 



(Sp(k, lu))/3 = r(k, u) 50ext(k, u) 

k 2 
— < 

47T 



e(k,u;] 



- 1 



^ext 



(k,w) 



defines the dielectric function e(k,u). Hence, 



k 2 

g(k,iu n ) = — 



i(k, iu n ) 



(A17) 



(A18) 
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It is worthwhile noting that, reverting to ordinary units and restoring Planck's constant 
h, uj n = 27in/hj3. Hence in the classical limit H — > 0, which is equivalent to the high 
temperature limit (3 — > 0, u n — > oo if n ^ and only the n = term contributes to the 
Fourier sum defining the correlation function. Thus in this limit, the correlation function 
becomes independent of the imaginary times, with 

1 



e(k,0) 



e ik -( r - r ') . (A19) 



The dielectric function is related to the irreducible polarization function by 



4:71 

e(k,io) = l + —U(k,u), (A20) 



which gives 



(i ■ \ n(fc,ia; ra ) 
<,(*,-„) = TTWj — y (A21) 

The "sum of ring graphs" approximation is obtained by approximating the exact polarization 
function which appears here by its lowest-order, one- loop form II^ ^ . This function is obtained 
by neglecting the denominator correction above and setting 

n (0) (r-r',r-r') = ((p(r, r)p(r', r')) + ^ , (A22) 

where the superscript on the expectation value indicates that it is computed for free fields. 
The charge density operator which appears here has the field expression 

p(r,r) = ^ /s ^(r,r)^(r,r), (A23) 

where the sum runs over all the particle species with charge e s (and over all spin components). 
In the free field approximation, 

(p(r,r)p(v',r')) + )^ = E/' ((^(r, ^](^ r0) + )^ ((^](r, r)^ s (^ r0) + )^ , (A24) 

since the total charge vanishes and with it the expectation value of a single charge operator. 
The field expectation values that appear here are obtained from the generic forms 
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V>(r, r) ^(r', r'))^ = / ^ [1 ± n(p)] e «p-("V*<p)(~ ') , (A25) 



and 



^{V', T') ^(r, r)) (0) = / iM n(p)e ^(r-r') e -^(p)(-r') ? (A26) 



'/3 J (2tt) 3 

where the upper + sign is for Bosons, the lower — sign is for Fermions 



E(p) = p 2 /2M (A27) 



is the kinetic energy of the particle with mass M, and 



"(P) = e-Mto-* T x ( A28 ) 



is the momentum distribution for free Bosons or Fermions. Using the results (|A25 ) and 
( |A26|) , we have 



U(0)( f r-r') - V P 2 ( _W_ P *k.(r-r') f ( d P) -\E s (p+k)-E s (p)](r-r') 

1 ' ] ^s s J (27r)3 e J (2vr)3 e 

n s (p)[l ±n s (p + k)] , r>r', 

(A29) 

n,(p + k)[l±n a (p)], r'>r. 
The Fourier transform of this expression yields 

1 } ^ s S J (2ir) 3 iuj n -[E s (p + k)-E s (p)] K } 

The long-distance behavior of the polarization is controlled by the Debye wave number 

/4 = 4vrn(0,0) . (A31) 

In the one-loop approximation, 



12 



It is easy to confirm the validity of the results ( |A25 ) and ( |A26| ) : They obey the free field 



equations of motion, the equal-time commutation or anticommutation relations, and they obey the 



cyclic boundary condition (A4) generalized to include the number operator and chemical potential 



in an effective Hamiltonian. These conditions fix the results uniquely. 
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n D — 



(dp) dn s (p) 



(2ir)* dE s (p) 

-^E/^/^(P)^.E/^(n s )r- (A32) 

In the dilute gas limit, the Bose and Fermi distributions can be replaced by the Maxwell 
distribution 

n(p) ~ e -^[ E (p)-^ , (A33) 

and in this limit 

Kl = ^PE/s(ns)f. (A34) 

The Debye wave number is essentially a classical quantity, and in the dilute approxima- 
tion the zero-frequency classical dielectric function is given by 

e(M) = l + ^|- (A35) 
Placing this in the classical limit QA19Q of the charge density correlation function gives 

^'^jS^l^r^- <a36) 

The corresponding potential correlation function is thus 

ww^^/^fA^- (A37) 

The partial fraction decomposition 

pp i J ~ ~L2~ h.2 i K 2 (A38) 

expresses the potential correlation function as the difference of the Fourier representations 
of the Coulomb and Yukawa potentials, and so 

(0(^0(0)^ ~ 1 [1 - exp{- KD |r - r'|}] . (A39) 
p\r — r'\ 

In the dilute gas limit, n(p) <C 1 and the Maxwell distribution ( |A33| ) applies. Thus, in 
this limit, and for r > r', the result ( |A29|) reduces to 
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[E s (p+k)-£ s (p)] (r-r' ) e -pE. (p) 



(A40) 



Since the energies that appear in the exponentials are quadratic in the momentum p, the 
momentum integral is Gaussian which can be evaluated by completing the square. The same 
considerations apply for the other time order r' > r, and one finds that for a dilute gas 



(A41) 



where 



C s (k,r-r') 



exp 



2M, 



r — r 



1-|t-t'|//3" 



(A42) 



APPENDIX B: DENSITY CALCULATED BY REACTION RATE METHOD 

The method employed in the text to compute the nuclear reaction rate in a plasma may 
also be used to compute the average density of one of the species of the particles in the 
plasma. This we shall do here so as to illustrate the method in a simpler context and also 
to obtain some results that are needed for the work of the text.R 

Using the techniques of the text, it is easy to see that the average density of a particle 
with chemical potential /i may be expressed as 

(n) p = e^(r' = 0, U{fl) ) p \v' = 0) , (Bl) 

where the overall matrix element is a single-particle matrix element of the particle in question 
and now 

U{0) = (exp j- ^Vg0(r(r),T)j) , (B2) 



13 A standard treatment of some of the results that we obtain here appears, for example, in Fetter 
and Walecka (1971), Sec. 30. 
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with q the charge of the particle and r(r) the particle's coordinate operator which undergoes 
free motion in the imaginary time r. The only approximation involved here is the require- 
ment that the particle density (n)p being computed is dilute. The other components in the 
plasma may be dense. 

To illustrate the method, we shall, however, compute the case in which all the particles 
in the plasma are dilute so that the exponential in U(/3) may be expanded with only the 
first non- vanishing correction retained, 

'"'lT ( lT , ((n(T(T).T)n(r(T).T))_ ; ; 



M/5))^l + y^^r'((0(r(r),r)0(r(r),r)) + ) . (B3) 



Here 



(0(r,r)0(r',r)) + ) = J ^ e Mr-*) G fa T _ r ') , ( B4 ) 



p J (2vr 

with, as is appropriate for the dilute case 



C (k r r') - 1 V p--n(r-r') (^V lK°>(fc,tt; w ) 



where the lowest-order polarization function n^ ) has been calculated in Appendix A. The 
work in the text [c.f. Eq.'s ( |3.43| ) and ( |3.44| )] tells us that, with M the mass of the particle 



in question, 



M^ 3/2 



(r' = 0,-^|r' = 0)= ^— J , (B6) 



and 

(r' = 0, ( e *W^)]) + |r' = 0) = (JL^J ^ C(k, r - r') 



(B7) 



Here 



C(k, r-r') = exp 1~\t - r'\ [l - |r - r'|//?] j (B8) 

is periodic in the imaginary times r and r' in the interval , (3 as is the function Gp(k, t — t'). 
Since two such periodic functions may be expressed as Fourier series, the double time integral 
reduces to a single integral, 
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r drdr'f(r - r')g{r - r') = (3 f drf(r)g(r) , 
Jo Jo 



(B9) 



and we may write our result as 

3/2 



(n) 



l i+ 4l! dT I^ c ^ G ^}- < mo » 



To evaluate the correction which appears here, we set 



C(k, r) = 1 + [C(k, r) - 1] , 



(Bll) 



and note that [C(k,r) — 1] vanishes as k 2 for small k. Thus the long-ranged contribution 
of the potential correlation function is reduced in this second term, it is no longer singular 
at small k, and the potential correlation function may be replaced by its zeroth-order value. 
On the other hand, the replacement of C(k,r) by unity just picks out the n = Fourier 
mode of the potential correlation function. Therefore, to within an accuracy of order e 2 , 



(dk) /4tt\ 4vrn(°)(A;,0) 



(2k) 



+ / dr 
Jo 



(2tt) 3 \k 2 J fc 2 + 47rn(°)(A;,0) 
(dk) 



(2nY 



C(k,r) 



%) 2 ^(k,r) 



(B12) 



Recalling that 47rII(0, 0) = k 2 d and noting that only potential singular behavior of the 
integration at small k can give a contribution that is not of order e 2 , we see that to this 
order we may write 

(dk) /4tt\ 47rn(°)(A;,0) f (dk) /4tT 



/ 



(2tt) 3 \k 2 J k 2 + 47rn(°) (k, 0) J (2tt) 3 \k 2 J k 2 + k 2 d 
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Here 



(dk) /47T\ k 2 d 



(2tt) 3 \k 2 J k 2 + K 2 D 



= K D , 



(B14) 



while 



/S(^) 2 [ n(0) ^o)-n(°)(o,o)] =-s / o °°rfA;[n( )(A;,o)-n( )(o,o) 



(2nf 



d 1 

dk k 



r/3 roo 1 A 

8 dr dk-i-U (0) (k,T), (B15) 
Jo Jo k dk 
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where we have used polar coordinates, integrated by parts, and re-expressed the zero fre- 
quency component of the polarization function as a integral over imaginary time. The results 
that we have just obtained may be restated in the form used in Section 3.2 of the text, 
rP . r (dk) . rP . r°° 1 d 



(B16) 



which we reiterate is accurate to terms including order e 2 . 

For the second set of terms in Eq. ( P12 ), we again integrate by parts as was done in 
Eq. (MTgl) to obtain 



o J (2tt) 3 



C(k,r) - 1 



47T 



k 2 J 



Vn(°)(fc, 



8 r^.r4{^N^) n(o) (^)]-^ n(o) (^)} 



(B17) 



lo Jo k ^dk*- i dk 

The second term in the integrand on the right-hand side of this equation just cancels the 
contribution in Eq. ( |B15| ). For the first term in the integrand, we recall the results ( A41 ) 
and ( |A42| ) from Appendix A which give the dilute gas form 



n (0) (^r) = £/s (ns) (o) exp 



A- 2 



2M S 



(B18) 



The exponential here is of the same form as that in the definition ( [B8| ) of C(k, r), and the 
two combine to form a single such exponential which involves the reduced mass 



1 1 1 



(B19) 



Thus, 



r°° „ 1 d 
u Jo kdk 



SI dr I <ll,--^ r [C( y k,T)U^\k,' 
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2M, 
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" V M« 



(B20) 



We write 
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n 2 DjS = Ane 2 s (n s ) p , (B21) 
which is the contribution of the plasma species s to the squared Debye wave number, and 

X. = g, (B22) 

which is the thermal wave length of a particle with the reduced mass M s . With these 
definitions, our evaluation reads 

r/3 roc I A _ 

8 dr dk--i- C(fc,r)n (0) (A;,r) 
Jo Jo k dk l 

Collecting the results finally yields 

<">« = e * (^l) 3 ' 2 1 1 + 1* 2 ' 1 " " T^A- l ) ' (B24) 

This result is accurate including terms of order e 4 . The first correction is the classical 
correction, and it is of order e 3 . The second correction, of order e 4 , is of a quantum nature. 
It involves the quantum wave length X s which, with ordinary units, is proportional to h. 

The light electron has a thermal wave length A e that is much larger than the thermal wave 
lengths A s of the nuclear particles in the plasma. Thus, the electron provides the largest 
quantum correction in the result (P24|) . The large thermal wave length of the electron, 
however, also implies that the chemical potential factor exp{/?/i e } for the electron is much 
larger than those of the other plasma particles since electrical neutrality requires that half 
the particles in the plasma be electrons but their number density is proportional to A~ 3 . 
Thus it may well not be a good approximation to treat the electron distribution as a dilute, 
Maxwell-Boltzmann distribution. On the other hand, the small electron mass implies that 
wave numbers k that are very small in comparison to \J M s / (3 dominate in the electron 
contribution to Eq. ( |B20| ), and so for their contribution we may write 

C(k,r)~l, (B25) 

up to a correction of relative order m e /M s . In this case, only the zero frequency mode of 
the electron polarization function appears, and the contribution of this function can easily 
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= "o/?E s 4,A. (B23) 



be evaluated since the electron Fermi-Dirac distribution may be written as a superposition 
of Maxwell-Boltzmann distributions, 

I oo 

„,„ „ 1 = V (_l)"+l e -^t^W-Me] (B26) 

e [3[E e (p)-n £ ] _i_ — ^ ^ V / 

Thus the electron contribution to Eq. ( |B2(J| ) may be expressed as an infinite sum of the 
Maxwell distributions already used in the previous evaluation, and we have 

rB poo 1 A 00 rnB poo 1 A 

8 dr d^n 0) (^)=8y(-ir +1 / dr dk^Tl^ >(k,T;n(3) 
Jo Jo kdk ^[ Jo Jo kdk 



z n=l V " L e 

1 00 
= _i e 2 me g(_ 1 )«+l 2e «^ 
^ n=l 
1 2 

= --e 2 m e — 5 , (B27) 

where the extra factor of 2 accounts for the two spin polarizations of the electron. Using 
this more accurate treatment for the electrons and noting that 

a = -j— (B28) 

e z m e 

is the electron Bohr radius, we find that 

(n) = e 0» (M ) 3/2 f 1 + -(3q 2 n D - ^£ V 4 s \ s - -fit. — ^ 1 . (B29) 

It should be recalled that we consider different spin states as separate species; the total 
number density of a spin 1/2 species is given by twice this result. 



APPENDIX C: CLASSICAL PLASMA 

As discussed in the text, the correction due to the ionic component of the plasma cannot 
be evaluated using a classical treatment for the plasma unless its thermal wavelength is much 
less than the turning radius of the reacting particles. In the text, we extracted this limit 
from our general result for a weakly coupled plasma. Here we shall show how this result is 
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obtained when one starts with a formulation that is often employed in the literature. This 
formulation neglects the motion of the center of mass and takes the reaction rate to be 
proportional to the single-particle, relative motion quantum-mechanical expectation value 

V V cl 







exp {-/3 [H r + ei 0(rM 2 /M) + e 2 </.(fM 1 /M)]} 







oc 



(3 



exp 



dr 



ei 0(r(r)M 2 /M) + e 2 0(r (r)Afi/M) 



(Cl) 



where we have passed to the interaction picture on the right-hand side of the equation [c.f. 
Eq. (|2.46| )1 and used the normalized expectation value defined in Eq. (|3.46| ). Here classical 
statistical mechanics is used to evaluate the thermal expectation value. This formula is often 
used in the literature to obtain an effective potential for the quantum-mechanical relative 
motion, and then the tunneling problem which describes the quantum expectation value is 
solved for this effective potential. We shall content ourselves by showing that this procedure 
gives our previous limit for the first plasma correction when the plasma is dilute. Using the 
result (|A39|) for the classical thermal expectation value of the potential correlation function 



and expanding out to the first correction, we get terms of the form 
1 rP 



1 + 2 J dTdT ' { e l{^ T ) M */ M ) <f>(Hr')M 2 /M))^ + 

- [ drdr' le 2 n . , , 2 — ; — — 
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1 - exp{-K D \r(r) - t(t')\M 2 /M} 
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(C2) 



Expanding in powers of the Debye wave number, which is equivalent to expanding in the 
small parameter K£>r max , and adding up all the terms produces 



Kd f drdr' 
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(C3) 



This is precisely the classical plasma limit that was discussed in the text following Eq. ( ^.71| ) 
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APPENDIX D: REAL TIME TROUBLES 



The work in the text made use of thermodynamic, "imaginary time" methods. Here 
we shall compare and contrast this method with the formulation which employs purely real 
time methods. The real time method might appear have the advantage of displaying the 
dynamics of the reaction process in terms of simple physical pictures, such as that proposed 
by Carrero, Schafer, and Koonin (1988). We shall show explicitly, however, that this is an 
illusion. Terms in the real time formulation that apparently have a straight forward physical 
interpretation may, in fact, be completely cancelled out by other terms. In particular, we 
shall show that the correction found by Carrero et al. has such a cancelation and hence does 
not exist. 

To relate the two formulations, we shall pass to an interaction picture. This is done by 
partitioning the complete Hamiltonian H of the total system into a part H that describes 
the dynamics of the background plasma and the reacting particles, but with no interactions 
between the reacting particles and the plasma, and the remaining part Hi that the describes 
the interactions of the reacting particles with the background plasma, 



H = H + H 1 . 




The interaction picture is obtained by writing 




(D2) 



where 




(D3) 



involves the time-ordered exponential with 



H^t) =e iHot H ie - iHot . 



(D4) 



The interaction picture time evolution operator is unitary, 
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U + (ti,t 2 y = U + (ti,h)- 1 , (D5) 

and obeys the group property, 

U+fa, t 3 )U+(t 3 , t 2 ) = U+fa, h) . (D6) 

The general expression for the reaction rate was derived in the text and presented in 
Eq. (|2.10 ), which we repeat here for convenience: 

»+oo 



oo 



dt l(dT)(rt(r,t)ic(py 



g 2 [ + °° dt f(dr) e iQt (4(t, t) ^(r, t) ^(r, t) V> 4 (r, t) ^1(0) 4(0) c- 2 ( ) r., ( ) ; , 



oo 



(D7) 

Using the interaction picture decomposition (|D2|) for the statistical operator exp{— /3H} 
which is the weight in the thermal expectation value {•••)$ and also using this decomposition 
for the time dependence of K^(r,i) produces 

dt (dr) (f/ + (-2/5,0)[/- 1 (t,0)/C t (r,t)f/ + (t,0)/C(0)) . (D8) 

-oo J P 

Here the superscript I indicates that the time dependence is now governed by Hq and also 
that the statistical ensemble is now described (except for the full normalizing partition 
function in the denominator) by exp{— (3Hq\. This is the "imaginary time" formulation 
that was essentially employed in the text. 

In the real time formulation, one computes the thermal average of the square of corrected 
matrix elements and thus arrives at 

T = J + °°dt y'(dr)([f/|(t,-oo)/C t (r,t)f/|(+oo,t)" 

[u + (+oo, 0)£(0)E/+(0, -oo)] Y p . (D9) 

To prove that this is indeed the same as the previous result ( |D8| ) , we make use of the group 
property ( |D6| ) and unitarity ( |D5| ) to write Eq. (P9|) as 



dt (dr) (f/|(0,-oo)f/ + (t,0)- 1 /C t (r,t)f/|(t,0)/C(0)f/ + (0,-oo)) . (D10) 
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Using the cyclic symmetry of the trace which defines the thermal average, we encounter 

U+(0, -ooJe-^EZ+^O, -oo) = e-P Ho U + (-ip, -oo - i(3)U^(0, -oo) . (Dll) 

It is implicit in the interaction picture representation that the interaction is adiabatically 
damped at large times (the "ze prescription"). Thus we may replace the complex infinite 
time limit — oo — i(5 by the real limit — oo and use the group property (|D6|) to conclude that 

f/ + (0,-oo)e~ /3 ^°[/- 1 (0,-oo) =e- pHo U + (-i(3,0). (D12) 



Using these results in Eq. (PIO ) reduces it to the imaginary time form (|D8| ) and proves its 



equivalence with the real time form 

To illustrate the differences in the real and imaginary time formulations, we shall examine 
a simple model in some detail using the methods of both forms. In this model, all the 
interactions of the final, produced particles are neglected (which may be done if the energy 
release Q is large). Thus, in this model, 

H, = f{dx)\p 1 <f> + p 2 <l>\, (D13) 

where p\ and p 2 are the charge density operators of the initial particles labeled by 1 and 2. 
To simplify the exposition, we shall also assume that one of the initial particles, say 2, is very 
massive. Thus, except for its ordering position, the charge density operator p 2 is independent 
of time. To further simplify the model, we shall neglect the Coulomb interactions between the 
initial particles 1,2 and compute only the correction involving a single interaction between 
each initial particle and the background plasma, the term proportional to which we 
denote by T^. Finally, in our model we shall take the limit in which the background plasma 
is treated classically. In this classical limit, the wave number vector k of a plasma correlation 
function is generally neglected relative to an initial particle momentum p since these enter 
in the combination p + ftk. As we shall see, this neglect of the plasma wave number may be 
done and gives the classical limit except when energy corrections are involved which lead to 
frequency shifts depending upon AE/h. 
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We first compute the eie2 term in the imaginary time formulation ( P8|) . As discussed 
in the text, the corrections associated with the time evolution operator U + (t,0) are negli- 
gible. We neglect these here and compute only those associated with the statistical factor 
U+(-iP,0), 

Pi(ri,r 1 )p2(r2,r 2 ))^/Ct(r,t)/C(0)\ / . (D14) 

Since the interaction between particles 1 and 2 is neglected, the charge density operators p\ 
and p 2 commute. Hence, since p 2 is time independent, the time-ordered product of the two 
charge density operators is independent of r 2 . In view of the discussion of Appendix A, the 
correction ( pi4| ) thus involves 



i 



e(k, 0) 



e ik-(n-r 2 ) 



= (0( ri )0(r 2 ))* , (D15) 

where in the second line we have written the classical limit. We stress that, without any 
further approximation, this correction entails the dielectric permittivity of the plasma e(k, u) 
at zero frequency, u> = 0. That is, the correction depends only upon the static properties 
of the plasma. Now, in the classical limit, the plasma wave number k is neglected, and the 
correction entails only the total charge operators of the initial particles, 

y"(dri)pi(n,r) = Q a , J{dr 2 )p 2 (r 2 ,r) = Q 2 , (D16) 

which are time independent. In the thermal expectation value in Eq. ( pi4|) , these operators 
just measure the charge (with the Boltzmann statistics for the dilute initial particles which 
we always take), and so we may make the replacements Q\ — > e\ and Q 2 —>■ e 2 to obtain 

r 12 = e ie2 (3 2 (0(0)0(0));' T , (D17) 

where To is the free gas reaction rate. The classical correlation function for a dilute plasma 
gives, according to Eq. ( |A39| ), 
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/3 (0(0)0(0))^ = ^, (D18) 
and thus we arrive at the basic Salpeter correction, 

r 12 = ei e 2 /3K D r . (D19) 
Let us note for future reference that, as discussed in the Introduction, 

^static pol = _ ei62KD (D20) 

is just the static polarization energy of the plasma induced when the two initial particles 
are on top of one another. Thus the Salpeter result may be written as 

r 12 = -/3u^ tic pol r . (D21) 

We turn now to compute the e\e 2 term from the real time form ( p9|) , which gives, on 
reading from right to left, 

r+oo f r+oo 



/+QO r r+oo r 

dt / (dr) / dt x dt 2 \ (dri)(dr 2 ) 
-co J J— CO J 

l - ((0(r 1 ,t 1 )0(r 2 ,t 2 )) + )^/Ct(r,t)(p 1 (r 1 ,t 1 )p 2 (r 2 ,t 2 )/C(O); 
+(<f>(r 2 , t 2 )0(r l5 h))^ ( (p 2 (r 2 , t 2 )tf(r, t)) _ (^(n, ^(O^ 
+ (0(r 1; t^fa, t 2 )) p ( ( Pl (n, t^/C^r, *)) _ (p 2 (r 2 , t 2 ))C(0)] 

(r 1 ,t 1 )0(r 2 ,t 2 ))_)^(p 1 (r 1 ,t 1 )p 2 (r 2 ,t 2 )/C t (r,t))_/C(O))^ [, (D22) 



P 
i 

P 
i 

P 

where the notation (• ■ •)_ denotes the anti-time-ordered product. To simplify the work, we 
immediately take the classical limit in which the momentum transfer Tik imparted by the 
plasma to the heavy particle charge density p 2 vanishes. In this limit, the integration over 
r 2 in Eq. ( P22|) gives the total charge Q 2 (t 2 ) which is time independent save for its ordering. 
With the Boltzmann statistics that we use for the reacting particles 1,2, Q 2 must appear to 
the left of the creation operator ijj 2 and to the right of the destruction operator ip 2 , and in 
these cases it is simply replaced by e 2 . This restriction of the ordering of Q 2 (t 2 ) limits the 
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time interval over which t 2 runs to — oo < £2 < or — oo < £2 < t. The interval from to £ 
described by the time evolution operator £/+(£, 0) will be neglected here just has it has been 
previously. With the Boltzmann statistics, a similar time restriction obtains for p(ri,£i), 
— 00 < £1 < 0. Thus, the real time result (|D22|) reduces to 
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dt / (dr) / dti / dt 2 / (dr x ) 

J J — OO J — 00 J 



(r 1; £00(0, t 2 )) + ) p e 2 (/Ct(r, £)/C(0) Pl ( ri , £ x ))^ 

+ (0(0, £ 2 )0(n, £ 1 )}^e 2 (/Ct(r, £)/C(0) Pl (n, £i))^ 
+ (0(r 1; £00(0, £2))^ (pi( ri , £ 1 )/C t (r, £)/C(0))^ 

-((0(n, £00(0, t 2 ))_) /3 e 2 (pi(r 1; tO^Cr, £)/C(0))^ | . (D23) 

The terms that appear in the successive lines on the right-hand side of this equation appear 
to have a clear physical description. The first line corresponds to the correction that results 
when the initial particles 1 and 2 have a Coulomb exchange which is modified by the plasma 
in a dynamical (frequency or time dependent) fashion before the fusion reaction takes place. 
The last line give the same correction for the complex conjugate amplitude that enters into 
the squared matrix element. The second and third lines do not describe such second-order 
interactions which take place between the initial particles before the fusion occurs. They 
instead describe the square of first-order amplitudes. They correspond to processes in which, 
if a complete set of final plasma states is introduced between lO and /C, the plasma is left 
in excited states. In one of the amplitudes, this excitation is caused by particle 1, in the 
other, by particle 2. The trouble with the real-time formulation is that there is extensive 
cancelation between these terms that might appear to have a distinct physical meaning. 
To make a first combination of terms, we write 

(r, £)0(r', 0) ± ) = ±\e{t - H) ( [0(r, £), 0(r', £')] ) p + ~ ({0(r, £), 0(r', £')})^ , (D24) 

where the curly brackets denote the anticommutator, and 

(0(r, £)0(r', £')) = \ ( [0(r, £), 0(r', £')] ) + \ ({0(r, £), 0(r', £')}) . (D25) 
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The terms involving the anticommutator pieces in Eq. ( |D23|) all cancel. As discussed in 
Appendix A, the commutator has the thermal expectation value 

( [«'■ '>■«'■ 0} >, = I H / % (F ) 2 «*' ' (D26) 

Appendix A also shows that c(k, u') appears as the spectral weight which defines the dielec- 
tric function e(k,u) by the dispersion relation 

1 , (D27) 



e(k,uj) k 2 J 2tt lu — to' — it 

a result that we shall soon make use of. 

Making use of Eq's. ( |A25|) and (|A26| ) in the Boltzmann statistics limit in which n(p) 



is neglected relative to 1, and neglecting the plasma wave number k in spatial momentum 
terms, but not in the energies as is appropriate to the limit in which the plasma in treated 
classically, the first two lines of the right-hand side of Eq.( p23| ) involve 

J (drje^ (V4(r, t)^) Pl {v u fc)), = ci / ^n 1 (p)e-^ r+iE ^ t e iAE ^ , (D28) 



and the second two lines entail 



J (d ri y k ' ri ( Pl (v 1 ,t 1 )4(r, t)M0)) p = e i / J^rMP + k) e -^ r+ ^V A ^ , (D29) 
in which 

AE = E 1 (p + k)-E 1 ( P ) = ^±^- 

lvn\ 

~^-^ = Vl -k, (D30) 

mi 

where the approximate equality applies to our classical limit. In view of these results and the 
decomposition ( |D24j) , we see that the first and last lines of Eq. ( p23j ) contain the combination 

-7ii(p)+m(p + k) = -e-^ + e"^ (p+k) ~-/3AEm(p), (D31) 

where, as always, we use Boltzmann statistics for the initial reacting particles, and the last 
approximation is appropriate for our classical plasma limit. The first and last lines on the 
right-hand side of Eq. (|D23|) are evaluated with the time integrals 
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/ dt\ I dt2e(h 

J—oo J — oo 



t \ e -iu/(ti-t 2 ) e iAEti 
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AE-u'- ie 



-uj' 



(D32) 



The integration of the denominators which appear here over u' with the weight c(k, uj') is 
just the dispersion relation ( P27| ), and we see that the first and last lines of Eq. ( |D23| ) 
produce (3 times the factor 



^ d 2 ynp0l (vi) 



(dk) 4vr 
(2tt) 3 F 



e(k, vi ■ k) 



1 



+ 



e(fc, 0) 



(D33) 



We schematically indicate the resulting correction to the reaction rate in terms of an effective 
average value of this dynamical plasma polarization by writing 



-pdyn pol 
1 12 



y-dyn pol^ 



vi))r . 



(D34) 



This is essentially the result obtained by Carrero, Schafer, and Koonin (1988) (to order 
e\e2). It corresponds to using the e\62 term of the classical, dynamical plasma polarization 
energy 



U(t) = y(dr)p(r,t)r°\r,t) 



(D35) 



evaluated for 



p(r, t) = ei<5(r - v x t) + e 2 5(r) , 



at t — when the two initial particles are on top of one another, 



Ug* p q1 (vi) = e 2 0f(0, 0) + ei ^ ol (0, 0) 



(D36) 



(D37) 



This result, however, is not the whole story. The dynamical aspects of this result, 
which corresponds to the second-order Coulomb exchange between the two initial reacting 
particles as modified by the plasma, are completely cancelled by the second and third lines 
of Eq. ( |D23| ), which correspond to squares of first-order amplitudes that leave the plasma 
in excited states. These terms have exactly the same structure as the dynamical terms that 
have just been computed except that the time integration (p32[) is now replaced by 
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dt\ I dt 2 e 



-iu/fa-ti) iAEti 



1 



AE 



1 



+ 



1 



AE — uj' — ie —uj 



(D38) 



This gives rise to a correction which is the same form as that in Eq. ( |D33| ) except that the 
terms in the curly braces there are replaced by 



1 



e(fc, vi ■ k) 



- 1 



+ 



e(k, 0) 



(D39) 



The addition of these terms cancels the dependence upon the frequency-dependent e(k, vi-k) 
and yields precisely the Salpeter result (|D21| ). The moral to the story is that simple physical 
pictures should be augmented with a correct formal basis to ensure that correct physical 
results are obtained. In addition to explaining the discrepancy between our results and 
those of Carrero et al. (1988), this cancelation may serve as a general caveat with respect to 
approaches that formulate the problem completely in terms of a plasma-modified two body 
interaction between the fusing particles. We should note, however, that, in one circumstance, 
the plasma excitation terms which give rise to correction involving Eq. ( |D39| ) vanish. The 
frequency - uj - dependence of e(k, uj) for a classical plasma arises from a plasma particle's 
velocity v s in the form v s • k. Thus, if the typical velocity of the reacting particle is much 
less than the typical velocity of the plasma particle, |vi| <C |v s |, then e(k, Vi • k) may be 
replaced by e(k, 0), and the plasma excitation correction involving Eq. ( |D39| ) vanishes. This 
is the case when the mass of the plasma particle M s is much less than that of the reacting 
particle, M s -C M u the electrons in the plasma. 

APPENDIX E: DYNAMICAL CORRECTION SIZE 



Here we shall obtain the order of magnitude of the leading "dynamical" correction given 
in Eq. ( |3.30| ). To do this, we take the dilute plasma limit to evaluate the real-time correlation 
function that appears in Eq. (|3.30|) as 

-4vrl 2 



(0(O)0(r,t'))/3 



W -ikT 

(2tt) 3 



k 2 



(El) 
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where Yl(°\k,—t') is the plasma polarization function given at the end of Appendix A in 
Eq's. ( |A41| ) and ( |A42| ), but continued to real time. With both r - r' > and t - 1' > [as 



is appropriate to the operator ordering in Eq. (]E1|) ], this continuation is given by r — r — > 
i(t — t'), and we have 

rt°\k,-t')^J2/ s {n s )fC s (k,-t'), (E2) 

in which 

C s (k,-t')=exp{i^t'j , (E3) 

where we have neglected the correction in the exponential involving t'//3 since in our appli- 
cation this is of the negligible order 1/ ' (3Q. The wave number integration over k in Eq. ([El]) 
may be performed by using the representation 

- = j o sdse~ sk \ (E4) 

interchanging integrals, and completing the square to obtain 

W0 W ,, fl), = E/, MS V5F f «h (. - ^) exp {- 4(s _ ;: /2Mj) } . (E5) 

The correction appearing in Eq. ( |3.30| ) involves the difference of this expression for r/0 



and r = 0. To estimate the size of this correction, it suffices to use the classical approximation 
for the operator r(t') ~ pt' /M. Noting that the produced nuclei move much faster that 
the initial fusing nuclei, we see that the largest effect of this motion is revealed by using 
v/2MQ as described in Eq. (^20|) , with now M a typical value of the mass of a produced 



nucleus. We change variables with s = ut'/2M s to exhibit the correction as 
±z/3( ei + e 2 )e r (\m [0(r(f),f) - 0(0)])^ 



o 



±i/3(ei + e 2 )e r J2/ s {n s )fV^ I dt' 

J 



f 



2M, 



o 



udu |( M -z)- 3 / 2 exp|-^-^-| -u- 3 A . (E6) 
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The u integration is well defined, and the exponent in it involves Qt' which at most is of 
order unity.0 The overall factor of l/y/2M s appears to imply that the light electrons in the 
plasma give the largest contribution. However, for the electrons, the mass ratio m e /M in the 
exponent in the u integral is very small, and to leading order in this ratio, the exponential 
can be replaced by unity. Thus, the leading contribution from the electrons in the plasma 
is independent of the nature of the particle involved in the reaction, except for the overall 
factor of the particle's charge e r . Since the ± sign in Eq. (|E6|) is positive for the initial 
fusing particles 1,2 [from the H^(t)] and negative for the final produced particles 3,4 [from 
the V(t)], and since charge is conserved, e\ + e 2 — e 3 — e 4 = 0, we find that the leading 
electron contribution vanishes. For the remaining contributions of the ions in the plasma, it 
suffices to replace the u integral in Eq. ( |E6| ) by a (complex) constant of order unity. Since 
t ~ 1/Q, the remaining t integral is of order 

f dt 1 Vt> ~ Q- 3/2 . (E7) 
Jo 

Remembering that the Debye wave number has the form k 2 d = 47re 2 f3(n) p, and that the 
plasma frequency of the ions u p has a size given by u 2 ~ k 2 d / f3M, we find that the leading 
correction is in fact of the order 

z/3( ei + e 2 )e r j\m [0(r(tV) - 0(O)])/3 ~ Pe x e 2 K D |^\/^} ■ ( E §) 

Here we have neglected the niceties of mass and charge ratio factors, denoting all of the 
nuclear masses by the typical mass M, and replacing (ei + e2)e r by t\t2- 



14 The factor which replaces Qt' for the initial particles is much smaller in view of Eq. (3.19) and 
the discussion below it. 
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